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We present a comprehensive study of the dynamical properties of the quantum Heisenberg an- 
tiferromagnet on a triangular lattice within the framework of spin-wave theory. The distinct fea- 
tures of spin-wave excitations in the triangular-lattice antiferromagnet are (i) finite lifetime at zero 
temperature due to spontaneous two-magnon decays, (ii) strong renormalization of magnon ener- 
gies Ek with respect to the harmonic result, and (iii) logarithmic singularities in the decay rate 
Fk. Quantum corrections to the magnon spectrum are obtained using both the on-shell and off- 
shell solutions of the Dyson equation with the lowest-order magnon self-energy. At low-energies 
magnon excitations remain well-defined albeit with the anomalous decay rate Fk oc A;^ at k — > and 
Fk oc |k — QafT'^^ at k — > Qaf- At high energies, magnons are heavily damped with the decay rate 
reaching (2Fk/£k) ~ 0.3 for the case S — 1/2. The on-shell solution shows logarithmic singularities 
in Fk with the concomitant jump-like discontinuities in Re[ek] along certain contours in the momen- 
tum space. Such singularities are even more prominent in the magnon spectral function A(k,w). 
Although the off-shell solution removes such log-singularities, the decay rates remain strongly en- 
hanced. We also discuss the role of higher-order corrections and show that such singularities may 
lead to complete disappearance of the spectrum in the vicinity of certain k-points. The kinematic 
conditions for two-magnon decays are analyzed for various generalizations of the triangular-lattice 
antiferromagnet as well as for the XX Z model on a kagome lattice. Our results suggest that decays 
and singularities in the spin-wave spectra must be ubiquitous in all these systems. In addition, we 
give a detailed introduction in the spin- wave formalism for noncoUinear Heisenberg antiferromagnets 
and calculate several quantities for the triangular-lattice model including the ground-state energy 
and the sublattice magnetization. 

PACS numbers: 75.10.Jm, 75.30.Ds, 78.70.Nx 



I. INTRODUCTION 

Heisenberg antiferromagnet (HAF) on a triangular lat- 
tice has been a focus of much attention as one of the basic 
model systems in which geometric frustration and low 
dimensionality are expected to yield new physical phe- 
nomena. Although available experimental realizations 
of the triangular-lattice antiferromagnei^ii^i^ii^i^ are de- 
scribed by such a model only approximately, either due 
to anisotropics or because of additional interactions, the 
ideal nearest-neighbor Heisenberg antiferromagnet on a 
triangular lattice given by 

7i = J^S,.S,, (1) 

remains the principal reference point. 

The semiclassical S ^ 1 triangular-lattice HAF or- 
ders in the so-called 120° structure, see Fig. [T] Histori- 
cally, it was anticipated that enhanced quantum fluctua- 
tions destroy the long-range antiferromagnetic order for 
the spin-1/2 model.'' However, calculations of quantum 
corrections within the spin-wave theory have suggested 
that the 120° magnetic structure remains stable even for 
S = l/2;^ii2iiii The early numerical results for small 



clusters were less conclusive, some supporting magnet- 
ically disordered stateiS and some confirming the spin- 
wave results. ^'^ Since the quantum Monte Carlo suffers 
from the infamous sign problem when applied to frus- 
trated models, it is not until the Green's function Monte 
Carlo worki^ that the magnetically ordered ground state 
of the spin-1/2 triangular-lattice HAF has been gener- 
ally agreed upon. More recent series-expansion-'^^. and 
density-matrix renormalization group (DMRG) studies^ 
have confirmed the stability of the 120° spin structure for 
the case oi S — 1/2 and yielded the value of ordered mo- 
ments (S) « 0.20, close to the previous result 

Gradually, it has been recognized that the truly dis- 
tinct physics of the quantum triangular-lattice antiferro- 
magnet concerns its excitation spectrum and the thermo- 
dynamic properties. The anomalous behavior of the lat- 
ter has been discovered earlier by the high-temperature 
series-expansion study. The temperature dependence 
of such quantities as entropy or susceptibility, exhibits 
significant differences between the triangular- and the 
square-lattice models: upon lowering temperature down 
to about J/2 the square-lattice antiferromagnet shows 
strong signs of ordering, while the triangular-lattice one 
does not. More recently, developments in the series- 
expansion method have allowed to calculate the excita- 
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FIG. 1: (Color online) The ordered 120° spin structure of the 
triangular-lattice HAF. 



tion spectra of the noncoUinear spin systems directlyj i"'^^ 
Numerical results for the magnon band of the spin-1/2 
triangular-lattice HAF deviate substantially from the lin- 
ear spin- wave theory (LSWT), with the overall band- 
narrowing by ~ 50%, flattening at the top of the spec- 
trum, and extra "roton-like" minima appearing at some 
special k-points. These results are so different from the 
well-known square-lattice case, where the spectrum is 
renormalized only modestly and almost uniformly up- 
ward, that Ref. [l^ suggested that at high energies the 
elementary excitations are not spin- waves but spinous. 
This hypothesis was questioned by the subsequent stud- 
ies which have looked into the role of magnon interac- 
tions within the framework of the spin-wave theoryJ^i^ 
It was shown that the first-order 1/5* correction strongly 
modifies the LSWT spectrum in an overall qualitative 
agreement with the series-expansion data. Further de- 
tailed comparison of the series-expansion and the spin- 
wave spectra has confirmed their qualitative similarities 
and outlined remaining differences 

Simply by the virtue of observing such significant dif- 
ferences in the triangular- and square-lattice spectra, 
these works have explained the contrast in the thermo- 
dynamic behavior of the two systems. Since the spin- 
wave bandwidth of the triangular-lattice HAF is reduced 
to ~ J and the other features such as the "roton- 
like" minimum are at the energies « J/2r^ the thermo- 
dynamic quantities must be dominated by these short- 
wavelength features down to much lower temperatures 
than in the square lattice case where the bandwidth is 
W~2J and the spectrum has a rather benign shape. 

In the recent Letter^° we have focused on another dis- 
tinct feature of the spin-wave spectrum in the triangular- 
lattice antiferromagnet: intrinsic damping due to sponta- 
neous decays. Note, that although the spectrum renor- 
malization in the triangular-lattice HAF is significant, 
the truly dramatic qualitative difference from the square- 
lattice case is the finite lifetime of the excitations at 
T ~ 0. Aside from yielding a substantial damping for 
excitations in the most of the Brillouin zone, the first 
1/5 correction produces the logarithmic singularity in 



the imaginary part of accompanied by the jumps in 
the real part along some contours in the k-space. We 
have analyzed the origin of such singularities and have 
related them to the topological transitions in the decay 
surfaces of magnons, which are due to the saddle-point 
van Hove singularities in the two-magnon continuum. 
Whether such logarithmic singularities are the artifact 
of the 1/S approximation of the theory or are the true 
features of the spectrum was only partially addressed in 
our worki^ For the decays into magnons that themselves 
acquire finite lifetimes, singularity will disappear in the 
higher 1/5' order and thus is not "real". However, if the 
singularity is due to decays into stable excitations, we 
concluded that the singularity remains essential and the 
magnon damping should remain very strong.^*' 

We would like to mention that the numerical study 
did not observe jumps in the real part of the spectrum, 
nor it reported the damping.^- This can be used as an 
argument against any "physical" singularities. However, 
by design, the series-expansion method finds the spec- 
trum that is purely real. Even with such restrictions, 
there where certain k-points for which the numerically 
obtained et had a convergence problem, shown by large 
error bars in Ref. [isl . 

Altogether, the qualitative questions remain: why 
there is such a substantial difference of the excitation 
spectrum in the triangular-lattice HAF from the more 
conventional square-lattice case? How generic are the 
decays, singularities, and anomalously large role of spin- 
wave interactions? What happens to the spectrum near 
the singular points? It is the purpose of the present work 
to address these questions. First, we would like to pro- 
vide a consistent overview of the spin- wave formalism as 
applied to the triangular-lattice model in order to demon- 
strate the origin of the drastic differences between the 
spectra of the collinear and noncoUinear ordered mag- 
nets. Second, we will elaborate on our previous findings 
on magnon decays and demonstrate that they must exist 
in a wide variety of frustrated spin models. Finally, we 
would like to extend our study beyond the on-shell ap- 
proximation of the previous works and clarify the fate of 
the singularities in the spectrum. 

The rest of the paper is organized as follows: 

• Sec. |TT] provides a qualitative discussion of the ori- 
gin of the spectrum anomalies in the triangular- 
lattice HAF and suggests a broader framework for 
the subsequent results. 

• Sec. [Hi] gives a detailed description of the spin- wave 
formalism for a noncoUinear HAF. Here, we also ad- 
dress the controversy over two different results for 
the 0(l/5^) correction to the sublattice magneti- 
zation of the triangular-lattice HAF, which exist in 
literatureJ^iii 

• The first-order 0{1/S) quantum corrections for 
the spin-wave spectrum are considered in Sec. IIVI 
Here we discuss the characteristic long-wavelength 
and short-wavelength features of both the spectrum 
renormalization and the decay rates. 
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• In Sec. |Vl we discuss in detail the kinematic de- 
cay conditions for a generic single-particle spectrum 
coupled to the two-particle continuum. The rela- 
tion between the singularities in the renormalized 
spectrum and the topological transitions in the de- 
cay surfaces is established in this Section. 

• Sec. |VI] is devoted to the discussion of the off- 
shell solution of the Dyson equation in the complex 
plane. We show that in the strong-coupling regime, 
the solution may cease to exist and the magnon pole 
may disappear in the vicinity of the logarithmic sin- 
gularity. The results of the numerical off-shell solu- 
tion of the Dyson equation for the triangular-lattice 
HAF are presented and it is shown that the damp- 
ing remains very substantial. In this Section we 
also present the results for the spectral function for 
several representative k-point and the quasiparticle 
residues across the Brillouin zone. We show that 
the log-singularities become even more prominent 
in the spectral functions, complicating the conven- 
tional analysis. 

• In Sec. IVIII we discuss several examples of other 
models in which decays and singularities are ubiqui- 
tous, such as the easy-plane XX Z and orthorhom- 
bically distorted J- J' models on a triangular lattice 
as well as the XX Z model on a kagome lattice. 

• Sec. lVIIll contains final conclusions and Appendices 
are used to provide more technical details. 



II. QUALITATIVE DISCUSSION 

A. Cubic anharmonicities 

In the case of quantum Heisenberg antiferromagnets 
on bipartite lattices in D > 2, the spin- wave theory 
agrees extremely well with the available numerical data, 
even for the "most quantum" spin-1/2 square-lattice 
modeli^'22ii2^'2i Therefore, not only one needs to under- 
stand the apparent anomalies found for the Heisenberg 
antiferromagnet on a triangular lattice, but also to ex- 
plain why such anomalies are not present in the prop- 
erties of the same model on bipartite lattices. In the 
following we give qualitative arguments that the unusual 
behavior stems from non-coUinearity of spins in the or- 
dered state, which is, in turn, induced by geometrical 
frustration. 

In quantum magnets with collinear spin con- 
figurations, the interaction between spin-wave ex- 
citations is described by quartic and higher-order 
anharmonicities i^^^i^ In contrast, a general quantum 
system with non-conserved number of particles is ex- 
pected to have anharmonicities of all orders beginning 
with cubic terms, which describe interaction between 
one- and two-particle states. The common examples 
are phonons in crystals^ and excitations in superfluid 
bosonic systemsjSS. In quantum spin-liquids, such as spin 



ladders and various dimer systems, triplet excitations 
may also have three-particle interaction terms.^^'"^^ 

For the ordered spin systems, cubic anharmonicities 
correspond to coupling of the transverse (one-magnon) 
and longitudinal (two-magnon) fluctuations, which would 
require the presence of mixing terms between S*^ and 
gx,y gpjj-^ components. Such terms are absent in the 
collinear Heisenberg magnets due to remaining C/(l) ro- 
tational symmetry about the direction of the magnetic 
order parameter. On the other hand, in the noncollinear 
antiferromagnets spin canting produces coupling of the 
transverse fluctuations in one sublattice to the longitu- 
dinal ones in the others. As a result, this yields cubic 
terms in the magnon-magnon interaction i^°i^^i^^i'^^'^"^i'^^ 

Noncollinear spin configurations in the Heisenberg an- 
tiferromagnets can be induced either by an external mag- 
netic field or by frustrating interactions. In the former 
case, spin canting and cubic anharmonicities are small in 
weak fields. At higher fields, cubic interactions dominate 
and lead to spontaneous magnon decays above threshold 
field ff* i36i37,38 jj^ frustrated magnets, spin canting and 
cubic terms are substantial already in zero field and thus 
play the key role in the spectrum renormalization. 

The above discussion creates a broader view on the or- 
dered quantum magnets and their spectra. Magnets with 
collinear spin structures have excitations that are intrin- 
sically weakly coupled. Not only the energy is minimized 
by the collinear spin orientation, but the interactions are 
also weak. Having a noncollinear spin structure necessar- 
ily implies much stronger coupling among the excitations. 
Thus, the collinear antiferromagnets should be consid- 
ered if not as an exception, but at least as a simplified 
subclass of quantum antiferromagnets. Naturally, many 
of their properties that are commonly assumed to be valid 
for all ordered antiferromagnets, such as the ubiquitously 
close agreement of the harmonic theory with numerical 
results, should not be expected to hold in general. 

Somewhat more formally, the lack of cubic anhar- 
monicities in collinear spin systems results in the absence 
of certain class of diagrams in the perturbative expan- 
sion. In the first order in 1/5, the only non- vanishing 
contribution is given by simple "balloon" Hartree-Fock- 
like terms (see Fig. |3]). Since these are w-independent, 
the corresponding correction to the spin-wave energy is a 
trivial redefinition of the interaction strength. The first 
tj-dependent correction appears only in the second order 
in l/S* and is already quite weak due to small phase-space 
factor. In noncollinear spin systems, cubic anharmonic- 
ity generates "bubble" diagrams already in the lowest 
1/5-order, not only inducing a substantial spectrum cor- 
rection but also providing a channel for the decays if the 
decay conditions are fulfilled. In the triangular-lattice 
HAF both effects are amplified because of the lower di- 
mensionality and because the tilting angle between spins 
is not small. Hence, the strong coupling of the longitu- 
dinal and transverse modes in noncollinear magnets in 
general and in the triangular-lattice HAF in particular is 
the key to understanding strong renormalization of their 
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spectra from the results of the harmonic approximation. 



of spin S the decay rate Fk is logarithmically enhanced 
relative to a perturbative result Fk ~ 1/5*. 



B. Decays and singularities 

Another important property of noncoUinear magnets 
is the ubiquitous propensity of the excitations to sponta- 
neous decays. The presence of transverse-to-longitudinal 
coupling is a necessary but not a sufficient condition for 
decays. In addition, the energy and momentum must 
be conserved in the decay process, yielding kinematic re- 
strictions. Decays depend, therefore, on the shape of the 
single-particle dispersion that may, or may not, allow 
for spontaneous decays. While we will give a detailed 
classification of various kinematic conditions later (see 
Sec.|V|, there are simple arguments for the decays to ex- 
ist in the triangular-lattice HAP. The complete breaking 
of the SO{3) rotational symmetry by the 120° spin struc- 
ture leads to three Goldstone modes fS, one at the center 
of the Brillouin zone and two at the ordering vectors ±Q. 
The former mode can be related to an infinitesimal twist 
of spins about the axis which is perpendicular to the spin 
plane, while the other two correspond to twists about the 
axes lying in the spin plane. The velocity of the latter 
modes is smaller than that of the former one. Such a 
difference guarantees that the energy conservation can 
always be satisfied for the decay of the fast quasiparticle 
into a pair of slower magnons, similar to the decay of the 
longitudinal phonon into a pair of the transverse ones in 
crystals.^* Since cubic anharmonicities necessarily gener- 
ate couplings between all magnon branches, this ensures 
finite lifetime for spin excitations in an extended part of 
the Brillouin zone at T = 0. The above consideration 
is trivially generalized to other frustrated antiferromag- 
nets where symmetry guarantees existence of more than 
one type of the Goldstone mode. We will elaborate on 
the decays in the XXZ anisotropic and orthorhombi- 
cally distorted triangular-lattice AFs as well as on the 
kagome AF in Sec. IVIII We would like to note that the 
kinematic conditions for decays can be completely sup- 
pressed by magnetic anisotropics, still leaving spins in a 
noncoUinear configuration. In that case decays are absent 
but the renormalization of the real part of the spectrum 
due to spin-wave interaction remains substantial. 

In the following, we also discuss extensively the singu- 
larities that occur in the decay rates of single-particle 
spectra due to the van Hove singularities in the two- 
particle continuum. While the main conclusion is that 
the essential singularities are cut-off either by the finite 
lifetime of the decay products or by the non-singular fi- 
nite lifetime of the decaying particle, the decay rates re- 
main parametrically enhanced by such singularities. A 
qualitative estimate of such a singularity-enhanced de- 
cay rate in the two-dimensional case is 



{v,ys)HAs/vi) 



(2) 



where V3 is the strength of a three-particle decay vertex 
and A is the momentum cut-off. Even for large values 



III. SPIN- WAVE FORMALISM 
A. Bosonic Hamiltonian 

The spin- wave theory of coUinear antiferromagnets on 
bipartite lattices is commonly formulated by introducing 
bosonic operators according to the number of magnetic 
sublattices.^-!^- In the Heisenberg triangular-lattice an- 
tiferromagnet spins form the three-sublattice 120° struc- 
ture at T = 0, see Fig. [TJ In order to go beyond the 
linear spin-wave analysisj^i^ it is convenient to trans- 
form to a rotating frame with the z-axis pointing along 
the local spin direction. Then, one needs to define only 
single species of boson operators with the wave- vectors 
that belong to the full paramagnetic Brillouin zone (BZ), 
Fig. [21 We would like to note that even in collinear anti- 
ferromagnets the above approach has substantial advan- 
tage when calculating higher-order spin- wave corrections. 
Furthermore, such a "single-sublattice" procedure has no 
alternatives when considering generic, incommensurate 
(spiral-like) antiferromagnetic ordering. 

As a first step, we assume that spins lie in the x~z 
plane and perform transformation from the laboratory 
frame (xq, zq) into the rotating frame (x, z): 



5f cos( 
5,f sin6 



S'i sin I 
Sfcosi 



(3) 



where 0^ = Q • and Q — (47r/3, 0) is the ordering wave- 
vector of the 120° spin structure. The spin Hamiltonian, 
Eq. ^ , in the new coordinate system takes the following 
form: 



n = J 



n 



sfs^ 



cos 



9,){S^S^ + SfSj 



\ ( QZ QX QX QZ\ 



(4) 



where (ij) denotes, as usual, summation over the nearest- 
neighbor bonds. 

The bosonization of Eq. ^ is performed via the her- 
mitian Holstein-Primakoff transformation: 



2S 



a]ai 



(5) 



{St 



Sf ± iSf) with subsequent expansion of square 



roots to the first order in alai/2S. Such an approxima- 
tion is sufficient to calculate the 0(l/5^) corrections to 
the ground-state energy and to the sublattice magnetiza- 
tion and to determine the 0{l/S) correction to the spin- 
wave dispersion. The resultant spin-wave Hamiltonian is 
given by: 



n 



sw 



Ho + n2 + n3 + n4 + o{s-^), 



(6) 
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FIG. 2: (Color online) Left panel: the Brillouin zone of a triangular lattice, lines are representative cuts. Central panel: 3D 
plot of the linear spin-wave energy tJk in the triangular lattice HAF. Right panel: intensity plot of ojk- Note different velocities 
of the Goldstone modes and different symmetries of the dispersion near F and K {K') points, and the saddle point at M {M'). 



where 7i„ denote terms of the nth power in the orig- 
inal (Holstein-PrimakofF) boson operators aj and a^. 
The first term in this expansion is the classical energy 
"Hq = —^JS'^N. Terms that are linear in and a] van- 
ish automatically because the 120° spin-structure corre- 
sponds to a minimum of the classical energy. Note, that 
the Hamiltonian ^ also yields a series in powers of S, 
with Ho = 0{S^), ii2 = 0{S^), Ha = 0{S^^^), and 
7^4 = 0(5°), respectively. The subsequent treatment 
of Ti-sw consists of diagonalizing the harmonic part Ti.2 
exactly and, then, treating H3 and H4 as perturbations. 



As a result, the linear spin- wave Hamiltonian takes the 
following form: 

no+n2 = -^JSiS + l)N + Y,£k(bibk + l) . (11) 



In the harmonic approximation, spin- waves are noninter- 
acting bosons with the energy 



Bl = 3 JS'wk 



where we define the dimensionless frequency 



(12) 



1. Linear spin-wave theory 

Noninteracting magnons are described by the linear 
spin-wave theory (LSWT) or the harmonic approxima- 
tion. After the Fourier transformation, the quadratic 
part of Tisw reads 



1 



k 



Ak = iJS{l + i7k) , 5k = ^JS-i^ , 
where 7k is a sum over the nearest-neighbor sites 



kx \/3, 

cos cos K, 



(7) 



(8) 



Diagonalization of TL2 is performed with the help of 
the canonical Bogolyubov transformation: 



Ok = "k&k + Wk^L 



under conditions — = 1 and 



2 I 2 
"k + ^'k 



Sk 



(9) 



(10) 



^k 



= V(l-7k)(l + 27k) 



(13) 



The shape of Wk is shown in Fig. O The harmonic spec- 
trum of the triangular-lattice HAF has several distinct 
features. The intensity map clearly shows that the ve- 
locities of the Goldstone modes at F (k = 0) and K, K' 
(k ~ ±Q) points are different. They are given by 



■iJS 



V3 



,(0) 



iJS- 



(14) 



Another notable feature is the clear three-fold symmetry 
of the modes near Q points. Instead of the usual convex 
and isotropic form, the spin-wave energy at small k — 
k — Q is nonanalytic with varying convexity: 

£k ~ vqkil — a^pk) , where a^p ^ cos3(/3 . (15) 

The overall shape of the dispersion is also more compli- 
cated than in the square-lattice antiferromagnet. At the 
M point, the center of the BZ edge, tJk has a saddle point 
with the energy roughly half of the bandwidth. Thus, al- 
ready with the harmonic spectrum the thermodynamic 
response at intermediate temperatures should be quite 
different from that in the square-lattice antiferromagnet. 

We introduce now several non-zero Hartree-Fock av- 
erages that contribute to the spin-wave corrections of 
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many static and dynamic quantities of the triangular- 
lattice HAF: 

n=(a|ai), m.^{a\aj), A = {aiaj), S={aj). (16) 

In the harmonic theory these are expressed as linear com- 
binations of the 2D integrals 



ci 



E 



(7k)' 

LOU 



(17) 





(c) 



(d) 



FIG. 3: The lowest-order vertices that yield 1/S-corrections 
to the spectrum and 1/5^ contributions to the static proper- 
ties. 



with I = 0,1,2, see Appendix \X\ for further details. In 
particular, the linear spin-wave correction to the stag- 
gered magnetization, (S) = S — 5S, is given by: 



JS-^n^; -(2co- 2-1- ci) = 0.2613032. (18) 



energy- and spin-dependence of the three-boson interac- 
tions explicit, we define dimensionless vertices related to 
the original ones by 



(21) 



Fi,2(l,2;3) = 3iJ^/— ri,2(l,2;3) 



2. Spin-wave interaction: eubic terms 

The cubic interaction terms Tis in Eq. ([5]) have no ana- 
log in coUinear magnets. They originate from the cou- 
pling of local and spin components in Eq. In 
terms of the original boson operators ([5]) the cubic inter- 
action is given by 



^3 J J -"^sinie j-Oi) ala^{a]+aj)-a]aj{al+a,) 
(y) 

(19) 

For the coUinear spin structures, sin(6'j — di) = Q and the 
cubic terms vanish identically. 

Performing consecutively the Fourier and Bogolyubov 
transformations from a|, ai to bl^, 6k operators in 
we obtain the interaction terms expressed via the "new" 
bosons: 



^3 = E[^ri(q,k-q;k)&^6[_q6k 

k,q 



(20) 



3! 



r2(q,-k-q,k) 6i6Vq4 + h.c 



Generally, terms linear in 6q and &q also appear after 
the above substitution. They represent quantum correc- 
tion to the pitch angle of the spin helix'^'^'^"-''-'^ or to the 
spin canting angle in an external magnetic field 
Such a correction vanishes in the triangular-lattice HAF 
because the ordering wave- vector Q corresponds to a sta- 
ble symmetry point in BZ. 

The first term in Eq. (j20p describes interaction between 
one- and two-magnon states and is symmetric under per- 
mutation of two outgoing momenta. We call it the "de- 
cay" term, although the decay processes may be only 
virtual. The second term in Eq. ([20]) corresponds to the 
spontaneous creation of three magnons and we refer to it 
as to the "source" term. The source vertex is symmet- 
ric under permutation of all three momenta. Making the 



with ri^2 given by 

fi(l,2;3) = 7i(mi + vi){u2U3 + V2V3) -f 72(1*2 + ^2) 

X (U1U3 -I- V1V3) - 73(1/3 + V3){uiV2 + V1U2), 
r2(l,2;3) = 71 (Wi +Vi){u2V3 +V2U3) +J2{U2 +V2) 

X {U1V3 + V1U3) + 73(7/3 + V3){uiV2 + V1U2), (22) 

where Ui, Vi are the Bogolyubov parameters (jlOp and the 
function 7k is defined as 



7k 



ifsi 
3 V 



sm — Z sm — cos — fc, 

2 2" 



(23) 



The diagrammatic representation of the decay and the 
source vertices is shown in Figs. [3l^a) and[3I^b), respec- 
tively. The above form of Fi 2 coincides with the expres- 
sions used by MiyakeiiSi^ 



3. Spin-wave interaction: quartic terms 

The last term in the spin-wave Hamiltonian ^ repre- 
sents the quartic terms: 



TI4 



t 3 / + + 

aja-j + - [a'-aiaiaj -\- a^ajajai 



/- t t 



h.c. 



(24) 



After the Bogolyubov transformation and normal order- 
ing of b operators, the interaction (|24p is replaced with 



H4 SEi + 5H2 



Ha 



(25) 



where the first two terms are the Hartree-Fock correc- 
tions to the ground-state energy and to the magnon self- 
energy, Fig. [31(c) and EJd), respectively. The magnon 
self-energy 6H2 contains both the diagonal and the off- 
diagonal terms, while for the square- lattice HAF the 
anomalous off-diagonal terms vanish. The final term H4 
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describes two-particle scattering processes. Similar to the 
case of the collinear antiferromagnetSf^l this latter term 
yields only higher-order 1/ S corrections compared to SE^ 
and 5H2 and is, therefore, neglected in the present work. 

To derive the explicit form of Eq. P5|) , it is technically 
more straightforward to apply the Hartree-Fock decou- 
pling (fT6|) to ((24|) first and use the Bogolyubov transfor- 
mation in 5H2 afterward. This is analogous to the treat- 
ment of quartic terms for collinear antiferromagnetsi^ 
With the details of this derivation delineated in Appendix 
El we simply list the end result for the correction to the 
ground-state energy: 



(co-t-ci-2c2-l)^ - -{Cl-C2f 



o 

and to the harmonic spin-wave Hamiltonian: 



(26) 



1 



J2'i'%b^-^Bl'\b^b.^ + bl^bl) , (27) 



where the first-order magnon energy correction and the 
anomalous self-energy terms are expressed as 



-(4) 



3J 



7^ (4co + ci~ 5c2 - 4) 



(28) 



-7k(2 - 2co + ci + C2) - 2co - 2ci + 4c2 + 2 



and Bi;> =-^(l-7k)(ci+2c27k) , (29) 



respectively. 



4- Effective Hamiltonian 

With the triangular-lattice model discussed in detail, 
we would like to outline the structure of the spin-wave 
Hamiltonian for a generic noncoUinear spin system. After 
the Holstein-Primakoff and Bogolyubov transformations 
with subsequent renormalization, when necessary, of the 
classical configuration, the spin-wave Hamiltonian takes 
the form of a polynomial of terms with increasing number 
of bosons, Ho + 'H2 + 'H3 + ■ ■ ■, and decreasing power of S. 
Therefore, the spin- wave expansion to order 0{S'^) will 
always result in an effective Hamiltonian: 

k 

+ J2[^K^^H-^^^ + F^,^blbl^_^bl + h. 

k,q 

The spin-wave Hamiltonian for a collinear antiferromag- 
net is the restricted form of (|30p with no three-boson 
terms present. The dynamic interaction among magnons 
in such a system occurs only in the next order due to 
the four-boson terms that are substantially weaker. In 
addition, for the commonly studied case of the nearest- 
neighbor square- and cubic-lattice antiferromagnets, the 




(a) 



(b) 





(c) 



(d) 



FIG. 4: The lowest order normal, (a) and (b), and anomalous, 
(c) and (d), magnon self-energies generated by cubic vertices. 



anomalous quadratic terms in Eq. ([30|) also vanish, leav- 
ing only a benign energy renormalization. Altogether, 
the role of magnon interactions in collinear spin states is 
significantly less important than for noncoUinear ones. 

The effective Hamiltonian (|30|) provides a basis for the 
systematic perturbative calculations of various properties 
of the triangular-lattice HAF. Parameters of the Hamil- 
tonian (|30[) for the triangular-lattice case were derived in 
previous subsections: 



ek = ek + 4'\ 



1 



B 



(4) 
k ' 



(31) 



T4,q = ^ri(q,k-q;k), Fk,q = ^r2(q, -k-q, k), 
where ek - 0{S^), 4'*^ and B^^^'> - 0{S°), and Fi 2 

0(^1/2). 

Using the standard diagrammatic technique for bosons 
at zero temperature, we define the bare magnon propa- 
gator as 



Go^(k,w) = w-Ek + iO 



(32) 



Then, the lowest-order diagrams contributing in the or- 
der 0{S'^) to the normal and anomalous self-energies are 
shown in Fig. [4] and in Figs.[3Kc), (d). The corresponding 
expressions for the normal self-energies are 



Sli^(k,c.) 



Sff(k,c.) 



|ri(q;k)|^ 



2 ^ UJ — Sq — £k-q + iO 



E 



|r2(q;k)p 



uj + e„+ £k+q - iO 



(33) 



(34) 



From now on we shall use, for brevity, only two momenta 
in the notations for the cubic vertices ([22|) . The anoma- 
lous self-energies are calculated as 



S^2^(k,C.) 



ly r2(q;k)Ft(q;-k) 

2 ^ UJ + eq + Ek+q - iO 



, (35) 



^l^.^(k,^)4E^^™^- (36) 



Ek- 



iO 
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Finally, the two frequency-independent contributions to 
the normal and anomalous self-energies are given directly 
by the Hartree-Fock terms: 



and to SE3 



Ef^^^(k)=£f and Ef/(k) = -B. 



^(4) 



(37) 



of Figs.[3jc) and[3l^d), respectively. The structure of the 
lowest-order contributions (|33P "- p7p remain valid, with 
necessary modifications of vertices, for an arbitrary non- 
coUinear antiferromagnet. 



B. Static properties 

Previous works on the spin-wave theory for the 
triangular-lattice HAF have calculated the second-order 
1/S corrections to the ground-state energy and to the 
sublattice magnetization. Two approaches have been 
employed for calculation of the latter: numerical ex- 
trapolation of the response to small staggered magnetic 
field^— and a direct diagrammatic expansioufii Surpris- 
ingly, they have produced two different results: AS2 = 
0.011/(25') versus 0.027/(25*), respectively. While the 
latter approach has to deal with more singular higher- 
dimensional integrals, the former one relies on numerical 
extrapolation in a small parameter. Below, we resolve 
the controversy over the value of the staggered magneti- 
zation in favor of the Miyake's result^^ by following the 
diagrammatic approach of Chubukov et al^ and point- 
ing out a delicate issue with numerical evaluation of can- 
celing singularities under integrals. For completeness we 
also briefly discuss the ground-state energy correction. 



1. The ground-state energy 

The first two terms in the 1 /S expansion of the ground- 
state energy of the triangular-lattice HAF are given by 
Eq. (dl]). Here, we calculate the next-order correction 
resulting from the magnon interactions. The contribu- 
tion from the quartic terms has already been obtained in 
the course of deriving the spin-wave Hamiltonian and is 
given by Eq. (pS)) . Another correction of the same order 
is generated by the source vertex. Fig. [3)3: 



|r2(q;k)|^ 



3! Ek 

k,q 



Ek-Hq 



(38) 



Combining all these contributions together one finds in 
the second order in 1/5: 



1 



h , (14 + 2/3) 
5 



(25)2 



(39) 



where /2 = (1 — X^k "^k) and the constants Ia and I3, are 
straightforwardly related to SE^: 



h = (co-Fci-2c2-l)2 - -(ci-C2)2 



-0.254293 (40) 



E 

k,q 



r2(q,k)2 



- tJk4 



0.13785(1) 



(41) 



The above four-dimensional integral has been calculated 
by two different methods: the Monte Carlo integration 
and the finite-size extrapolation of lattice sums using 
clusters with different aspect ratios^^ both agreeing in 
all significant digits. Altogether, the ground-state energy 
in the second-order of 1/5 expansion is 



0.436824 0.02141 



25 



(42) 



The above result agrees with the previous calculation^ 
improving on the numerical accuracy of the last term. 



2. Sublattice magnetization 

To calculate the staggered magnetization we use the 
diagrammatic approach which is very close, aside from a 
few technical details, to the one used in Ref. [ij. Within 
the spin-wave approach the sublattice magnetization is 

(S) = S - (ala,) = S -SS , (43) 
where the quantum correction SS is expressed as 

= E H + ("k + «k) (^k^k) + 2ukWk(&k&-k)] . (44) 

k 

The first term under the sum is the LSWT result^ al- 
ready given by Eq. ^ 



0.2613032 . 



(45) 



The two remaining terms in (|44p contain bosonic averages 
which vanish in the LSW approximation and contribute 
only to the next order in 1/5. Therefore, we write: 



(5) = 5 - SSi - 



SS2 



25 ' 

where the last correction has two contributions: 

= (552,1 + <552,2 , 



(46) 



(47) 



6S: 



, CJk ^ , t^k 

k 



Calculation of the bosonic averages in the above ex- 
pression must be performed to the first order in 1/5. As 
explained in Appendix [Bl these averages are straightfor- 
wardly related to the normal and anomalous self-energies, 
see Figs. [3] and 21 In particular, the magnon occupation 
number (fek^k) is only due to Y]^^^ from Fig. Ub), while 



9 



the other two normal self-energy corrections, E^"-* and 
Sf^-^, have zero contributions. On the other hand, all 

three off-diagonal self-energies, S^2\ ^i?! ^'^'^ ^12^: con- 
tribute to (6k^-k)- Leaving the details of the derivation 
to Appendix iBl we present here the final answer: 



5S2 = 



9 , ^ 7k(l - 7k) 

k ^ 



16 



9y 7k y ri(q;k)r2(q;-k) 



(48) 



3v-- 1 



E 



(tJq + O^k-l-q + Wk)^ 



This expression agrees with the formula derived previ- 
ously in Ref. [ll| apart from the corrected sign in front of 
the third term. 

As is often the case with the higher-order spin-wave 
corrections, the individual contributions in Eq. ((^5]) are 
divergent: the integrands in the second and the third 
terms behave as 0{l/k^) at k — > Q, which means that 
not only the leading divergences in them, but also the 
sub-leading ones 0{l/k^) must cancel in order to pro- 
duce finite result. Expanding in small Sk = |k — Q|, 
such a cancellation can be verified analytically.^^ Still, 
the expression given in Eq. (|48|) is not well-behaved nu- 
merically. If one tries to evaluate 6S2 directly using the 
Monte Carlo integration, the outcome appears to be di- 
vergent. If some other methods are employed, the result 
may seem to be regular. We have used the simple finite 
sums in the k-space that correspond to periodic clusters 
with the subsequent finite-size extrapolationpi^ For any 
given subset of (rectangular) clusters with the fixed as- 
pect ratio the result of Eq. converges to a finite value 
as the size of the cluster L ^ 00. However, as an indi- 
cation of the problem, subsets with different aspect ratio 
yield different values of 682 in the thermodynamic limit. 

The origin of the problem is the following. The inter- 
nal integrals in Eq. ([^5]) over q are not divergent and, 
generally, scale with the lattice size asM^: 



E 



-Fk,q = * 



(00) 



L 



(49) 



In the thermodynamic limit, two such terms cancel near 
certain points and regularize the ~ 1/k^ singularity in 
the external integral over k: 



1 r 
fc3 



4- 



(00) 



(00) 



1 



[Ak^ 



(50) 



However, numerically such a cancellation is not complete 
as it carries a 1/L term as in ([^5]) : 



1 



(L) 



1 

fc3 



Ak' 



a 
L 



(51) 



Since the 2D integral of 1/k^ diverges as L, the 1/L cor- 
rection from mutually canceling terms in (j5ip will give an 



unphysical contribution to the L — > 00 limit of Eq. (|48p . 
This explains the erratic behavior of the numerical values 
of 6S2 and suggest that the extra care should be taken 
with Eq. The way to regularize this problem is 

described in Appendix |B] After the regularization, the 
final answer for Eq. (|48|) can be obtained by any standard 
integration method, which yields: 



SS2 = -0.011045(5) 



(52) 



This result differs from SS2 ~ 0.027 quoted by Chubukov 
et al.^^ , which may have suffered from the above integra- 
tion problem. Our value agrees, though, with 6S2 ~ 
0.011 obtained by Miyake^^ by a different method that 
avoids highly singular integrals but deals with extrapo- 
lation in a small auxiliary field. 
Finally, 



(S) - 0.261303 - 



0.011045(5) 
25^ 



(53) 



Thus, the spin-1/2 Heisenberg antiferromagnet on a tri- 
angular lattice has the following value of the ordered mo- 
ments: (S) w 0.24974 within the second-order spin-wave 
expansion. Note, that this spin-wave value is somewhat 
larger than the results of the Green's function Monte 
Carlo,— the series-expansion;^ and the recent DMRG 
calculation&lS that give (S) « 0.205(15). 



IV. SPIN- WAVE SPECTRUM: 1/S 
CORRECTION 

Similarly to the calculation of the static properties, 
perturbative expansion for the magnon spectrum has to 
be performed order by order in 1/5 to ensure cancella- 
tion of all possible divergences in the individual contribu- 
tions. This difficulty notwithstanding, derivation of the 
ffi'st 1/S correction to £k is straightforward. The anoma- 
lous terms do not contribute in that order and the new 
pole of the magnon Green's function is determined by 



e - ek - E^^ (k) - (k, s) - S^/ (k, e) = , (54) 

to which we refer to as to the Dyson equation. The self- 
energies in (I54p are given by the diagrams in Figs, 
and [Ha,b) . Solving it self-consistently in the complex 
plane for a new renormalized spectrum £ = Ck — iT^ con- 
stitutes the off-shell approximation discussed in Sec. IVII 
The first-order 1/S correction to the spectrum is ob- 
tained within the so-called on-shell approximation. In 
this approximation the self-energies are evaluated at the 
bare magnon energy e = £k- This leads to the following 
expression for the renormalized spectrum: 



(b), 



£k - «rk = £k 



q 



-(4) 



ri(q;k)^ 



(55) 



r2(q;k)= 



LJk- 



Wk — *0 



" ^k+q + ^k 
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In the right-hand side the harmonic energy is 0(5*), while 
the rest of the terms are 0(5"). The finite magnon decay 
rate Fk comes only from the first term in the brackets. 
Since the Goldstone modes should be well-defined in the 
ordered AFs, we expect Fk ^ Ck in the long-wavelength 
limit. The details for that limit are given in the following 
subsection, while Sec. IV. B is devoted to the behavior of 
the renormalized spectrum in the full BZ. 



A. Low-energy magnons 

1. Velocity renormalization 

The triangular lattice antiferromagnet has three Gold- 
stone modes at k = and ±Q. The existence 
of these zero-energy modes follows directly from the 
broken 50(3) rotational symmetry in the spin space 
and, therefore, should not be affected by quantum 
renormalizationsi^S As a consistency check of the 1/S ex- 
pansion, it is important to verify the presence of acoustic 
modes in the renormalized spectrum (I55p . Such a ver- 
ification was first performed in Ref. 11 where the 1/5 
corrections to the velocities of the Goldstone modes were 
derived. We have reproduced these corrections with an 
improved numerical accuracy, although, technically, a 
somewhat different route was followed for the derivation. 
Numerical values of the spin-wave velocities are: 



VQ 



3JS 



3J5 




V3 / 0.11488 

0.08285(2) 
25 



(56) 



2. Long-wavelength decays 

In the long-wavelength limit decay rates can be always 
calculated perturbatively because of the smallness of in- 
teraction among low-energy excitations and due to re- 
duction of the phase-space volume available for decay 
processes. In other terms, presence of the well-defined 
Goldstone modes implies smallness of the damping rate 
with respect to their energy, Fk ^ Ek- As with the ve- 
locity renormalization, such a behavior has to be verified 
using Eq. ([55]) . As was remarked before, existence of the 
cubic interactions alone does not immediately yield fi- 
nite lifetime for the excitations since the corresponding 
decays may be virtual. The decay processes become real 
if both the total energy and the total momentum can be 
conserved simultaneously. Since the two-particle excita- 
tions form a continuum of states, the energy conserva- 
tion can be rephrased as the requirement of an overlap 
of the single-particle branch with the two-particle con- 
tinuum. This yields certain kinematic conditions on the 
bare spectrum Sk, which are discussed in Section V us- 
ing the example of the triangular-lattice HAF. Here we 
assume that such conditions are fulfilled and consider 
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FIG. 5: (Color online) Spin- wave decay rates in the vicin- 
ity of r and K points in the triangular-lattice HAF. Dots 
are numerical solutions of Eq. (|55p and dashed lines are the 
asymptotic results (see text). 



several scenarios for the decays of the long-wavelength 
excitations relevant to the present model. 
Generally, the decay rate is given by: 



E 



(57) 



where Vk_q is the decay vertex. Due to the energy con- 
servation, the upper limit for the momentum of created 
quasiparticles should be of the order of k. Then, for the 
linear spectrum £k ^ fc, a- naive power counting suggests 
the following answer 



„D-1 



(58) 



where comes from the D-dimensional phase space, 
1/fc is due to reduction of that space to decay surface 
from the energy conservation, and 14 is the typical am- 
plitude of the cubic vertex on the decay surface. Let us 
assume that the cubic interaction follows the standard 
form^^ Vk.q oc \/kqq', q' = |k — q|, and that in a typi- 
cal decay process the final momenta are q,q' ^ k. This 
makes |Vfep cx k^ and yields a seemingly universal power 
law for the decay rate: 



Fk 



k^+^ 



(59) 



For D — 3 this yields Fk ~ k^ which matches the result 
for the decays of the convex phonon branch, but, as we 
shall see shortly, only coincidentally. 

In reality, the situation is more delicate and possible 
power-law exponents for the decay rate asymptote de- 
pend on the specifics of the problem. In the case of a 
single weakly nonlinear acoustic branch, relevant to the 
phonon spectrum in ^He, 



eic = ck + ak 



(60) 



the decays are allowed only for a positive curvature of the 
spectrum, a > (convex Sk). In this case, the unstable 
quasiparticle emits two excitations in a narrow solid angle 
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FIG. 6: (Color online) Spin- wave energy for the S = 1/2 triangular-lattice HAF along the symmetry directions in the BZ. The 
dashed and solid lines are the results obtained in the LSWT approximation (sk) and with the first-order 1/S correction (eic), 
respectively. The vertical arrows indicate singularities and intersection points with the two-magnon continuum. Gray areas 
show the width of the spectral peaks due to damping Pk from Fig. [T] Dotted lines represent the minimum of the two-magnon 
continuum obtained from the LSWT spectrum. 



centered in the direction of the initial momentum. The 
apex angle of the decay cone scales a.s 6 ^ k such that the 
phase space is k^^~^ instead of . Then the restriction 
from the energy conservation gives l/kO"^ ^ 1/k^ instead 
of 1/fc in the previous consideration. Altogether, for the 
case of the cubic upward curvature of the spectrum, the 
answer is universal (see also Ref. |42): 



HAFs in ZJ-dimensions: 



(61) 



which yields Fk ~ k"' in 3D29'44 ^nd Fk ~ k^ in 2D. The 
2D result applies to the square-lattice HAF in a strong 
magnetic field where the convexity of the sound-like 
branch changes from a < below the threshold field H* 
to a > in the high- field region H* < H < if^'^'^'^ 

If several acoustic modes with different velocities are 
present, the fast quasiparticle can always decay into two 
slow ones. This situation is simpler than the previous 
case since the nonlinearity of the spectra plays no role. 
The phase space factor becomes k^~^ now, in agreement 
with the above naive consideration. Physical realizations 
of this scenario include decays of the longitudinal phonon 
into two transverse ones in solids^^ as well as the decay of 
the k ^ into two k ^ ±Q magnons in the triangular- 
lattice HAF. Clearly, such a channel of decays withstands 
quantum renormalization of the velocities and is perti- 
nent to the other noncollinear AFs. 

Interaction between phonons in crystals obeys the con- 
ventional scaling asymptote Vk,q oc ^/kqq' and, conse- 
quently, |Vfcp cx k^. Therefore, the naive power counting 



of Eq. ([59|1 is valid for them. However, in the case of the 
triangular-lattice HAF, the result is yet different from 
([55)1 because the three-magnon vertex (PO)) is anomalous 
and scales as Vk,Q+q oc {q' — q)y/k/qq' for small k. For 
a typical decay process q,q' ^ k giving: |Vk,Q+qP ~ k 
instead of k^ . Altogether this yields for the noncollinear 



■D 



(62) 



Direct analytic expansion of Eq. ([55)1 gives the following 
decay rate of the k ^ Q magnons in the triangular-lattice 
HAF {D ^2): 



Ek 



9PJ 
4\/27r 



TT 

V32 



1 



2\/2 



(63) 



which is Fk ~ 0.789 Jfc^, in agreement with the numerical 
data in Fig. [H 

The kinematic conditions for magnons near ±Q [K 
and K') points are more subtle. The magnon energy 
has a nonanalytic expansion in small k = k — Q with 
varying convexity: £k ~ c/c -|- a^pk'^ , where the nonlin- 
earity depends on the azimuthal angle ^ cos3(/9, see 
Eq. (jlSp . This form, together with the commensura- 
bility of the AF ordering vector Q with the reciprocal 
lattice of the crystal, allow for momentum and energy 
conservation for magnon decays from the steeper side 
of the energy cone at k ^ Q onto the flatter sides at 
q, q' — > — Q. Commensurability of Q is important here 
as 3Q = is the necessary condition for the conservation 
of the quasi-momentum. Thus, magnons near the Q- 
point are unstable only in a certain range of angles. The 
decay vertex for k — > Q magnon has the conventional 
scaling: VQ+k,-Q+q oc \/kqq', so the decay probabihty is 
~ k^. However, due to a constraint on the angle 
between k and q, the decay surface in q-space is a cigar- 
shaped ellipse with length ~ k and width ^ fc'^/^. That 
makes the restricted phase volume of decays to scale as 
^(3D-5)/2^ This results in a nontrivial law: 



Fk ^ 



(64) 
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FIG. 7: (Color online) The spin-wave damping of the 5 = 1/2 triangular-lattice HAF calculated in the first-order of the 
1/S expansion along the same lines as in Fig. [G] Arrows indicate singularities and intersection points with the two-magnon 
continuum. 



which gives k''/'^ for the decay rate in 2D. Numerically, 
along the KT hne Tk ~ 1.2Jk'^/^, see Fig. O Away 
from this direction the damping exhibits an anomalous 
angular-dependence Ft ~ l/(cos3(p)'^/^. Such a behav- 
ior is related to the saddle-point singularities which is 
discussed in the next Section. 

B. High-energy magnons 

The renormalized spin-wave energy and the magnon 
damping given by Eq. ((55|) have been calculated using 
the Monte Carlo integration method with 10^-10^ inte- 
gration points in the full BZ. Numerical results for the 
spin-1/2 case are shown in Figs. [S] and [7] for represen- 
tative symmetry directions in the BZ, see Fig. [2] for no- 
tations. Statistical errors of the calculation are compa- 
rable/smaller than the corresponding line width. The 
first prominent feature of the spectrum is that the renor- 
malization is stronger at large momenta, in a qualita- 
tive agreement with the series-expansion resultsJ^i^^i^ 
As we have argued before, this is a consequence of the 
considerable coupling between single-magnon excitations 
and two-magnon continuum determined by cubic anhar- 
monicities. The momentum dependence of the minimum 
of the continuum is shown in Fig. [5] by dotted curves. 
The intersection points of these curves with the magnon 
branch are marked as fc^'s. Inside the continuum, magnon 
excitations become intrinsically damped acquiring a non- 
zero Fk, see Fig. [71 The magnitude of the overlap of the 
continuum with the single-particle branch gives a quali- 
tative idea of the phase space available for decays. The 
damping rate is also illustrated in Fig. [S] as the shaded 
area between Sk — Tk and £k + Tk- 

Another interesting property of the renormalized spec- 
trum in Fig. [S] is its shape in the vicinity of the M-point 
(edge center of the BZ). Quantum renormalization con- 
verts the saddle point of at kjv/ into a local minimum 



surrounded by flat parts j ' i Such local extrema must 
contribute significantly to the thermodynamic properties 
of the triangular-lattice HAFJ^ii^ The minimum in ek 
is more pronounced in the numerical results^^ than in 
Fig. [51 It was called a "roton-like" minimum and it was 
suggested that it might be a signature of spinonsj^S This 
has been questioned later since the same feature in the 
spin-wave results can be explained by the enhanced den- 
sity of two-magnon states near the M-point 

Although the above discussion is important for bench- 
marking the spin-wave theory with the series expansion 
results, the most remarkable property of the renormalized 
spectrum is the substantial jump-like discontinuities in 
Ek, marked as k* points in Fig. [51 with the jump heights 
reaching 1/4 of the magnon bandwidth. If considered 
without the concomitant behavior of Fk, such jumps in 
Ek are especially enigmatic. The values of the damping 
at the top of the magnon band are also quite substan- 
tial, leading to broadening of the spectral peaks with the 
widths about (2Fk/ek) ~ 1/3. The most striking fea- 
tures of Fk are the sharp logarithmic singularities at sev- 
eral points along the selected cuts of the Brillouin zone. 
These are precisely the same k* points where the jumps 
occur in the real part of the spectrum. Thus, the two 
singularities are intrinsically related. Analytic consider- 
ation of these singularities in the magnon self-energy is 
presented in Sec. V.B. Here we simply state that their ori- 
gin is due to the intersection of the single-magnon branch 
with the line of the van Hove saddle-point singularities 
in the two-magnon continuum. 

The logarithmic singularities in the damping rate sig- 
nify a breakdown of the 1 / S'-expansion in the vicinity 
of the singular fc*-points. It is remarkable that the log- 
singularities are obtained already in the lowest-order in- 
teracting spin-wave theory. Therefore, if any other prop- 
erty calculated with the 1/ S'-expansion is to be compared 
with numerical methods, such singularities must be un- 
derstood and their proper treatment within the theory 
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FIG. 8: (Color online) The decay contours in q-space for 
magnons with selected k's along the TK direction; k — k' and 
k = k* are the same as in Fig.[7l left panel. The corresponding 
q-contours undergo topological transitions at k' and k* . 

outlined. Note, that the numerical series-expansion re- 
sults exhibit lack of convergence (large error-bars) at cer- 
tain k-pointsii^ Although, this is not the same as singular 
jumps, such numerical features might occur due to sim- 
ilar reasons. Analytic results for the spectrum near the 
singularities and the proper treatment of the latter will 
be presented in Sees. IVl and IVI[ respectively. 

Altogether, the main effects of quantum fluctuations 
in the spectrum of the triangular-lattice HAF are the 
substantial magnon damping in the major part of the 
BZ, singular behavior of the decay rate along certain 
contours in the k-space, and strong downward renor- 
malization of the magnon energies compared to the har- 
monic spin-wave theory. All these effects underline the 
importance of cubic anharmonicities in the noncoUinear 
AFs and represent major qualitative differences from the 
collinear cases. 



stable with Fk = 0. The region in the k-space with sta- 
ble excitations is separated from the decay region by the 
decay threshold boundary. Generally, two-particle exci- 
tations form a continuum of states in a certain energy 
interval 

< i?k(q) = eq + ek-q < Er"" ■ (66) 

Thus, the decay threshold boundary is determined by 
the intersection of the single-particle branch with the 
bottom of the continuum E^^'^. For our problem, the 
decay region is the hexagram shown in Fig. [5] 

Needless to say, the two-particle decays may be pro- 
hibited in the entire BZ, a situation that is common for 
collinear antiferromagnets in zero field.— In such a case, 
spontaneous n-particle decays with n > 2 are also prohib- 
ited since all the energies in the n-particle generalization 
of Eq. ([65|) are positively defined. Thus, the analysis of 
the two-particle decay conditions is crucial for determin- 
ing whether any spontaneous decays are possible. In the 
case when the two-particle decay region is finite, the n- 
particle decays can, in principle, be allowed in a wider 
region of the k-space, see, e.g., the corresponding study 
for the phonon branch in ^Hc;^ However, such an expan- 
sion of the decay region requires rather special conditions 
and is not discussed in this work. 

Aside from finding whether spontaneous decays ex- 
ist or not, there is another important reason for con- 
sidering the decay kinematics. As was emphasized by 
Pitaevskii22ii^ in the analysis of the two-roton decay 
threshold in superfluid *He, the enhanced density of two- 
particle states near the minimum of the continuum may 
produce strong singularities in the single-particle spec- 
trum at the decay boundary. This yields various unusual 
effectSf^i^i^ including complete disappearance of the 
quasiparticle branch inside the continuum. Apart from 
the singularities in the spectrum at the decay boundary, 
additional singularities may occur within the decay re- 
gion due to topological transitions of the decay surfacci^ 
General analysis of these two effects and its application 
to the triangular HAF are discussed here. 



KINEMATICS OF TWO-PARTICLE DECAYS 



A. Decay threshold boundary 



The aim of this section is to consider kinematic con- 
straints that follow from the energy conservation in the 
two-particle decay process: 



(65) 



This equation should be treated as an equation in q with 
the initial momentum k as a parameter. The solutions of 
Eq. (|65p form the decay surface in the q-space. Examples 
of the decay surfaces (contours in 2D) for the triangular- 
lattice HAF are shown in Fig. [5] for a few representative 
k-points along the TK-line. 

As a function of k, the decay surface changes and may 
disappear completely. In such a case, particles become 



The boundary of the decay region corresponds to the 
intersection of the single-particle branch with the bot- 
tom of the two-particle continuum ([66]) . To find one 
should analyze the extrema of the continuum given by 



5Sk(q) _ de. 



dq 



q 



de. 



q 

dq 







(67) 



k-q 



Thus, the extremum condition is satisfied if the velocities 
of the decay products are equal. This yields D equations 
with q and k as independent variables. To find the decay 
boundary one needs to solve the decay and the extremum 
conditions and (IST)) together. Solving them yields a 
set of {D — l)-dimensional surfaces in the k-space. Since 
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FIG. 9: (Color online) The Brillouin zone of a triangular lat- 
tice. The shaded area shows the region where spontaneous 
two-magnon decays are allowed. The lines correspond to the 
extrema in the two-magnon continuum described in text. 

the non-trivial solution may exist not only for the mini- 
mum of the continuum but for the other extrema as well, 
some parts of these surfaces define the decay threshold 
boundary while the rest of the solutions will correspond 
to the special surfaces within the decay region. The latter 
are considered in Sec. fVlB. 

For the gapless spectrum, one should also verify 
whether the emission of the acoustic (Goldstone) exci- 
tation at q ^ Qi corresponds to This condition 
is separate from Eq. (|67|) as the Goldstone modes are 
not the extrema but the end-points of the spectrum. In 
the case when one or more of these gapless modes define 
-Ejf' the decay boundary is given by 

£k = £k+Q, • (68) 

Although the above discussion covers two most gen- 
eral cases, there exist a particular form of the solution 
of ((65)) and (|67p for which finding the threshold bound- 
ary simplifies greatly. The two-particle continuum al- 
ways possesses an extremum when both decay products 
have equal momenta, up to a reciprocal lattice vector G: 
q, k — q = (k ± G)/2. The condition of equal veloci- 
ties (1571) is automatically satisfied in this case. Such an 
extremum crosses with the single-particle branch on the 
surface determined by 

Ek = 2£(k+G)/2 • (69) 

Among the previously studied cases, the decay bound- 
aries for the square- and cubic-lattice HAFs in a strong 
external field satisfy the above equation. 

For the case of the triangular-lattice HAF, let us begin 
with the analysis of the decay involving the Goldstone 
modes. The acoustic branch emerging from q = does 
not correspond to the crossing point of £k and E^"^ and 
should be disregarded. The condition ((68)) on the emis- 
sion of Goldstone modes with q = ±Q can be rewritten, 
using the expression (|12p for the harmonic spectrum, as 
7k = 7k^Q- This last equation is easily solved and the 



results are shown in Fig. [9] by the dotted lines. One can 
show that q = ±Q points correspond to the absolute 
minimum of the continuum for all k within the shaded 
area in Fig. [5] To demonstrate that, ek±Q is also plotted 
in Fig. [6] as dotted line. 

We begin the analysis of the extrema in the two- 
particle continuum with the decay threshold involving 
magnons with equal momenta. The solutions of Eq. ()69p 
for the reciprocal lattice vectors G1.2 = (±27r, 27r/\/3) 
and G3 = (0,47r/\/3) are readily found numerically and 
are shown in Fig. [9] by solid lines. No solution exists for 
G = 0. Solid lines in Fig. [S] lie entirely within the shaded 
area and correspond to the saddle points of the contin- 
uum, not the minima. The solution for the general situ- 
ation of the threshold decay into a pair of non-identical 
magnons with equal velocities but different momenta and 
energies, ek-q 7^ Eq, is obtained by solving numerically 
the decay ((65p and extremum ((67p conditions simultane- 
ously. The result is plotted in Fig. (H) by the dashed line. 
As in the previous case, this contour corresponds to the 
line of saddle points of the continuum within the decay 
region. 

Altogether, the decay region is given by the union of 
all areas obtained in the above three cases. As one can 
see, the decay threshold boundary in the triangular HAF 
is determined entirely by the emission of acoustic e±Q 
magnon. In accordance with the long- wavelength consid- 
eration of Sec. IIVI A. the area around k = is precisely 
where the decay of the fast magnon ek^o into two slow 
ones Eq^iQ takes place and this area is completely en- 
closed in the decay region. Similarly, there are only finite 
segments in the vicinity of the ±Q-points where decays 
are possible. The other thresholds due to two-magnon 
extrema do not contribute to the decay boundary. As 
is shown in Sec. IVIIl such a mutual arrangement of dif- 
ferent thresholds may or may not be the case for other 
closely related systems. Already a simple addition of 
the XX Z anisotropy modifies the decay boundary and 
switches saddle points into minima for certain momenta. 

The damping of magnons close to the boundary of the 
decay region can be considered in a manner similar to 
the long-wavelength analysis of Sec. IIVI A. 2. Since the 
velocity of the decaying magnon near the boundary must 
be larger than the velocity of the emitted acoustic ±Q 
mode, the phase space factor is (k — k\j)^~^ as is given 
in Eq. ((55)) . The probability of the decay has a smallness 
only with respect to g ~ (fc — kb). This yields the 2D 
decay rate near the threshold due to acoustic mode: 

Fk - (fc - hf (70) 

in agreement with the results in Fig. (7) in the vicinity of 
kb points. 

B. Topological transition in the decay surface 

A truly remarkable feature of the magnon spectrum 
in the triangular-lattice HAF is the logarithmic peaks 
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in Fk obtained in the first-order of 1/S expansion. Af- 
ter preceding discussion it is clear that the location of 
such singularities corresponds to the crossing of single- 
magnon branch with the surface of saddle-points of two- 
magnon continuum. Indeed, the A;*-points in Figs. [6] and 
[7] belong to the threshold contour for the decays into 
pairs of equivalent magnons, solid lines in Fig. [5] Much 
weaker anomalies, which are visible only as small peaks 
in Fk that are denoted by A;'s in Fig. [71 correspond to 
the threshold contour for the decays into non-equivalent 
magnons, dashed lines in Fig. [31 

In the vicinity of the crossing points of the single- 
particle branch with the saddle-point surface of the con- 
tinuum, the magnon decay surface undergoes a topologi- 
cal transition, see Fig.[8l Expanding the energy conserva- 
tion condition ([H5)) in small Ak — k— k* and Aq — q— q*, 
where k* is the point on the threshold contour and q* is 
the saddle point of _Ek* (q) , we obtain 

Ao2 Ag2 

Ek - eq - Ek-q ~ (Vl - V2) • Ak - ^ + -J^ = , (71) 

where Vi and V2 are the velocities of the initial and final 
magnons, and a, h are constants. Depending on the sign 
of (vi — W2)AA: the solution of the above equation are 
two conjugate hyperbolas, which transform into a pair of 
crossing lines for Ak = 0. Thus, k = k* corresponds to 
the point where the decay surface splits into two disjoint 
pieces, as is indeed observed in Fig. [51 

The discussion in terms of saddle points of the con- 
tinuum and topological transitions of the decay surface 
complement each other. While the van Hove singular- 
ities are always present in the two-magnon continuum 
(|66p they do not necessarily cross with the one-magnon 
branch. On the other hand, using the topological tran- 
sition perspective one can argue that the occurrence of 
such a crossing does not depend on a precise form of 
Ek- As was discussed in Sec. lIVl A. acoustic magnon with 
small k decays into two magnons that are close to Q 
and — Q. Its decay surface consists of two disjoint parts 
near K and K' points, see = tt/A in Fig. [51 Simi- 
larly, a magnon with k ^ Q can emit two magnons in 
the vicinity of — Q. The corresponding decay surface is 
a single closed contour near the K' point, see = tt in 
Fig. [51 Therefore, moving along an arbitrary trajectory 
in k-space between F and K points the decay surface 
must undergo at least one topological transformation, 
which also implies crossing of at least one saddle-point 
threshold contour. The actual number of such transfor- 
mations is determined by the short-wavelength details of 
the spectrum. For the triangular-lattice HAF there are 
two of them as is demonstrated in Figs. [51 and [51 

Let us now consider behavior of the magnon self-energy 
in the vicinity of singular points. The decay vertex is 
regular at k — > k*, q ^ q* and gives an unimportant 
constant factor. Using the expansion (|7T|) we obtain for 
the singular part of the magnon self-energy: 



A straightforward integration in (|72p yields 



E(k,ek)« J 



ReS(k, Ek) ~ sign(AA:) , Fk ~ In 



A 

aF 



(73) 



[Fk = — ImE(k, £k)]- The cut-off parameter A is deter- 
mined by the characteristic size of the region in the k- 
space where the expansion (|7T|) holds. The linear size of 
the smallest "droplet" of the decay surface at the topo- 
logical transition can be taken as an upper bound on 
A, see Fig. [8l Such an estimate explains the difference in 
the strength of anomalies in Fk for k* and k' points. The 
topological transition at k* consists of joining/splitting 
of the two approximately equal contours of substantial 
size, while the k' point corresponds to the splitting off of 
a small piece. 

To put this discussion in a broader context we note 
than in the earlier works^2i^d^ the situation was consid- 
ered when a singularity occurs at the boundary of the de- 
cay region rather than in the interior. In such a case, the 
extremum in the two-particle continuum that is crossed 
by the single-particle branch is a minimum, not a saddle- 
point, and the analog of Eq. ([72]) for I](k, ek) is given 
by 



E(k,ek)cx j 



{vi ~ V2)Ak - ql/a"^ - q'^/b'^ + iO' 



(74) 



After integration this yields the following characteristic 
anomaly 



ReE(k,ek) ^ In 



A 

Afcl 



F. 



e(Afc) 



(75) 



{vi - V2)Ak - ql/a'2 + q^/b^ + iO' 



(72) 



where 9(x) is the Heaviside step function. Thus, the 
situation is reversed in comparison to our case: the log- 
anomaly occurs in the real and the jump in the imagi- 
nary parts of the spectrum. Since the imaginary part of 
I](k, £k) is related to the two-particle density of states, 
in 2D it is natural to have a jump in Fk upon entering 
the continuum and a log-singularity upon crossing the 
saddle-point line inside the continuum. By the Kramers- 
Kronig relations such jumps and logs in ImE(k, ek) re- 
sult in logs and jumps in Rcl](k, ek), respectively. For 
the 3D systems, logarithmic peaks disappear and one ob- 
tains only square-root singularity: ReE(k, Sk) — \/\Ak\. 

Another important question concerns whether singu- 
larities in the spectrum will survive the higher-order 1/5* 
treatment. If the singularity persists, vertex corrections 
may become important, see Appendix [Cl^^ Our Sec. IVll 
discusses this problem. 



VI. OFF-SHELL DYSON EQUATION AND 
SPECTRAL FUNCTION 



The unusual logarithmic singularities in the magnon 
decay rate Fk found in the first-order 1/S corrections 
signify a breakdown of the standard spin-wave expan- 
sion. They represent an extra theoretical challenge and 
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FIG. 10: (Color online) Renormalization of the decay region for the spin-1/2 triangular-lattice HAF. Left and middle panels: 
Dashed and dotted lines are the LSWT predictions for the magnon dispersion £k (dashed) and for the minimum of the two- 
magnon continuum (dotted); solid and dotted-dashed lines represent same results for the spectrum Ek with the first-order 
quantum correction, respectively. Right panel: Modified decay region; lightly shaded regions highlight the area of uncertainty 
due to the re-entrant behavior seen in the middle panel. 



have to be renormalized in order to obtain the actual 
dynamical behavior. The purpose of this Section is to 
describe, at a technical level, different approaches to this 
problem. In Sec. IVII A we show that the singularities 
are regularized, for the most part, by allowing for the 
finite lifetime of the magnon in the initial state within 
the so-called off-shell approach. However, in the strong- 
coupling case, the single-particle excitation may disap- 
pear in the vicinity of the singularity, similar to the ter- 
mination point in the quasiparticle spectrum of super- 
fluid *He.'*^ In Sec. IVII B we discuss the magnon spectral 
function ^(k, lj) and find additional singularities in the 
w-space that are directly connected to the van Hove sin- 
gularities in the two-magnon continuum. 



A. Singularities and off-shell solution 

1. Modified decay region 

The first-order quantum correction to the magnon dis- 
persion (Sec. IIV[) leads to a significant narrowing of the 
magnon bandwidth for the spin-1/2 triangular-lattice 
HAF. Therefore, one may ask how strongly the spectrum 
renormalization affects the decay boundary and whether 
the decay condition Eq. (|65p is still satisfied for the renor- 
malized £k. In the long- wavelength limit, the difference 
between the velocities of the acoustic modes vq > vq, 
see Eq. ([55]) . guaranties that the kinematic conditions 
for decays at k — > are always fulfilled. While this ar- 
gument is applicable to any system with several types 
of the Goldstone mode, the decay boundary is defined 
by the short-wavelength features, specific to a particular 
system. 

In Fig. [TO] we present numerical results, which show 
that the size of the decay region does not change ap- 
preciably even for the S = 1/2 triangular-lattice HAF. 



The left and the middle panels demonstrate the renor- 
malization of the bottom of the two-magnon contin- 
uum for two representative paths in the BZ. The lowest- 
energy two-magnon states within the decay region still 
correspond to the emission of the ±Q Goldstone mode, 
^min _ e^^Q. While the renormalizations of the spec- 
trum and the continuum are significant, the intersection 
points of the two change only weakly. For the YM path 
in the BZ one observes a re-entrant crossing of the single- 
magnon branch with the bottom of the continuum, which 
leaves some uncertainty in defining the new decay bound- 
ary. This behavior is related to the jump-like singulari- 
ties in the single-magnon e^, which should become more 
well-behaved once the singularities are regularized. The 
right panel of Fig. [TO] presents the "new" magnon decay 
boundary for the spectrum that includes l/S renormal- 
ization. The light-shaded regions show the uncertainty 
areas where the re-entrant behavior of the spectrum and 
the continuum occur. Overall, the decay region does 
not change significantly in comparison with the LSWT 
boundary in Fig. [5] 

2. Higher-order diagrams perspective 

One way to regularize the decay diagram is to allow 
for a finite lifetime of the decay products by dressing the 
inner lines in the "bubbles" in Fig. H] The effect of such a 
dressing depends, however, on whether the saddle-point 
momenta of the decay products q* and k* — q* in Fig. [5] 
fall inside or outside of the decay region. If q* lies inside 
the decay region for at least one of the final magnons, 
then this magnon will acquire a finite lifetime in the next 
order. The logarithmic singularity in Fk will be removed 
in this case since the energy conservation law Eq. (|65p is 
now satisfied only on average. For the triangular-lattice 
HAF such a scenario is realized for a large fraction of the 
"weak" singularities {k' points in Figs. [HJ [7]). However, 
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all of the "strong" singularities {k* points in Figs. [H [T]) 
and some of the "weak" ones belong to another class in 
which the saddle points for both magnons created in the 
decay process are outside of the decay region. Hence, 
at the saddle points, the logarithmic divergence of the 
one-loop diagrams will persist even for the renormalized 
spectrum and the singularities seems to remain essential. 

Thus, the above approach requires summation of an 
infinite series of diagrams that contain leading-order di- 
vergences, similar to the Pitaevskii's treatment of the 
spectrum termination problemj ^^i^^ Such a treatment is 
hindered within the spin- wave theory by the divergence of 
the individual terms at k — > Q in each order of the l/S 
expansion, the problem already mentioned in Sees. IIIII 
and IIVI While a qualitative statement on the result of 
such a regularization can be made (see Appendix [C| any 
quantitatively reliable calculation are problematic in the 
light of this problem. 



3. Off-shell Dyson equation 

The above seemingly hopeless situation is resolved if we 
note that the logarithmic singularity occurs for a magnon 
which is inside the decay region. This means that the 
Dyson equation ([M]) should be solved "off-shell," i.e. 
with the same complex energy e — e^ — iT^ both outside 
and inside E(k, e). Physically, such a procedure allows 
for a finite lifetime of the initial magnon while magnons 
created during the decay process remain stable. The off- 
shell approach avoids complications related to multi-loop 
diagrams and, as we demonstrate below, is sufficient to 
regularize the singularity. Note that, the magnon energy 
in the off-shell solution contains corrections beyond the 
first 1/5' order. Nevertheless, it can be shown that it is 
yet free from the k ^ Q divergences associated with the 
higher-order diagrams mentioned above. 

After a methodological remark made in AppendixlPjon 
the proper sign of the imaginary part of e in the decay- 
like self-energy, the Dyson equation is given by 



£-ek-S(k,e*) = , 



(76) 



where e* is the complex conjugate of e and S(k, e) 
includes all one-loop contributions given by Eq. (|54p . 
Rewriting the above equation explicitly for the real and 
imaginary parts one obtains the following system 



ek = ek + Re[S(k,ek + irk)] , 
Pk = -Im[E(k,ek + «rk)] >0 



(77) 



Let us first demonstrate how a finite Pk can regularize 
the singularity. Replacing in Eq. I|7ip with complex 
e and using parametrization s — Sc — ■ Ak = p — 
|p|e~*''', where Sc is the position of the saddle point in 
the continuum, we obtain after integration 



Here, oc rf(k*,q*) J'^S and A is the momen- 
tum cutoff. Omitting the non-singular contributions, 
Eqs. (|77| are rewritten as 



Ek — £k 



JS 



AJS 



Ek 



V2 • Ak -I- iPk 



(79) 



The on-shell solution is recovered by substituting £k = 
Ec + 'Vi ■ Ak instead of e = £k + «Pk in ([79]) . This yields 
p = (vi — V2) • Ak and also implies that cos ip = ±1 
(or 93 = 0, tt) depending on the sign of Afc. As a result, 
one finds the jump in the real part of the spectrum and 
the log-singularity in the decay rate in agreement with 
Eq. dZSl). 

Any renormalization should shift the "bare" singu- 
larity into a new crossing point of the single-magnon 
branch with the surface of the saddle points in the two- 
magnon continuum. Assuming that the real part of the 
energy renormalization is already included in the defini- 
tion of Ec and k* and that Ek can be still expanded as 
Ek — Ec + vi • Ak we obtain 



AJS* 



Av • Ak -f iTu 



(80) 



As a result, the imaginary part of the solution at the 
singular point Ak is now regular and can be determined 
by solving the transcendental equation: 
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r 



(81) 



with Pk = Vij* /JS. Depending on the relative strength 
of the three-magnon coupling {V3/ JS)"^ to the size of the 
dimensionless momentum cutoff A there are two different 
regimes: 



Pk 
Pk 



Vi 



In 



JS 
AJS 



A{JS)^ 



^3 



Vi/A{JSf « 1 , 
VilA{JSf > 1 . (82) 



Thus, at large V3^/A(J5')^ the decay rate is independent 
of the coupling and is defined by the phase volume factor. 
Recalling that Vi cx J^S" in the case of the triangular- 
lattice HAF, we obtain estimates: 

^ociln(.A) 



Pk 

JS 



cx A , 



5A> 1 , 
S'A < 1 , (83) 



where the first expression is relevant to the "strong" sin- 
gularities with large phase space volume for decays (A:*- 
points), and the second is for the "weak" ones (/c'-points). 
For k ^ k* , one finds for the off-shell solution 



Vi 





+ z In 


AJS" 














P 





(78) 



cos Lp 



Re(Ae - V2 • Ak) ^ Av Ak 
lAe- V2 • Akl Pk 



. (84) 
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FIG. 11: (Color online) Comparison of the magnon spectrum (upper row) and the decay rates (lower row) obtained in the 
harmonic approximation (dotted lines), with the first-order 1/S corrections (dashed lines), and by solving the Dyson equation 
(solid lines), respectively. 



Hence, — s- tt/2 and the singular jump in et, see 
Eq. ((79|l . also disappears in agreement with the above 
assumption. 

In a hypothetical case of the strong cubic term 
(V3 S> JS) one needs to consider vertex renormaliza- 
tions discussed in Appendix [C] Briefly, this replaces 
the log-singularity in the self-energy with ImE(k, e) cx 
l/ln|A/p|. Solving the Dyson equation yield the same 
answer as in ([82|) for the /lV{JSY < 1 limit, while in 
the opposite case solution for the single-particle spectrum 
near k* does not exist. This is similar to the complete 
disappearance of the spectrum at the termination point 
of the spectrum in the superfluid ^Hc;^. 

Finally, we present in Fig. [Tl] the numerical solution 
of the Dyson equations (|77|) for the spin- 1/2 triangular- 
lattice HAF. While the jumps and the logarithmic 
peaks disappear, the damping rate remains substantial 
throughout BZ. Note, that the overall shape of the off- 
shell £k is in a better agreement with the series expansion 
data^^ than the on-shell results. In particular, the "ro- 
ton" minimum is well-pronounced and the "flat regions" 
are much less significant than in the latter case. How- 
ever, there is an overall upward energy scale offset of 
our results relative to the numerical ones. This may be 
due to both the remaining higher-order l/5'-corrections 
to the spectrum in the spin-wave theory approach and 
the neglect of the imaginary part of the spectrum in the 
series-expansion calculations. The upward energy renor- 
malization of the off-shell versus on-shell results is natu- 
ral as the latter tends to overestimate the energy shifts. 

Altogether, the main result of the off-shell considera- 



tion is that this approach naturally resolves the singu- 
larity problem, regularizing the log-singularities in the 
decay rates and removing the concomitant jump-like dis- 
continuities in the real part of the spectrum. As the 
result, the decay rates remain significant and are loga- 
rithmically enhanced relative to the perturbative results. 

B. Spectral function 

Generally, a detailed information about both the pole- 
like and the incoherent parts of the single-particle spec- 
trum is obtained from the diagonal component of the 
spectral function 

A{\i,Lo) ^ --lm\Gii{\i,uj)\ . (85) 

TT 

In the leading one-loop approximation the diagonal 
Green's function for the triangular-lattice HAF is given 
by 

Gri'(k,c^) =Lu-ek-S(k,u;) , (86) 

with i;(k,w) = S^^(k) + X;li'(k,w) + Sfi^(k,w) ex- 
pressed by Eq. (|54p . In quantum antiferromagnets the 
spectral function A(k, w) is also related to the dynam- 
ical structure factor 5'(k, w) which is directly measured 
in inelastic neutron experiments. Generally, 5'(k, u) also 
has contributions from the off-diagonal and two-particle 
correlations^MSiill but the spectral function ([55)1 still 
provides the major component. 
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In the absence of intrinsic damping, the quasiparti- 
cle peak in ^(k, w) occurs precisely at w = ek found 
from the solution of the Dyson equation. In the pres- 
ence of spontaneous decays the solution of Eq. ([77)1 dif- 
fers from the position and the width of a quasiparticle 
peak in the spectral function because the latter is de- 
fined on the real w-axis. Another characteristic feature of 
A(k, Lo) for all non-coUinear AFs is the contribution from 
the two-magnon continuum due to a non-orthogonality 
of the one- and two-particle excitations. In particular, 
for the momenta k inside the decay region, the spectral 
weight in A(k, uj) should become non-zero above the bot- 
tom of the continuum. On the same ground, one should 
also expect singular behavior due to the van Hove singu- 
larities of the continuum to be prominent in the spectral 
function at any k, not only at special contours of k*. 
This is because frequency scans through all possible en- 
ergies and is not limited to the "mass-surface" lo = e^. 
Since we are restricted to the one-loop approximation for 
S(k, uj) due to difficulties with the higher-order diagrams 
discussed above, such van Hove singularities will appear 
as sharp features in A(k,w). This is due to both the 
one-loop approximation for the self-energy and because 
u! at which the system is probed is purely real. Qualita- 
tively, all singularities are expected to be regularized by 
the higher-order contributions. The complication is, of 
course, that a quantitative calculation of such a regular- 
ization can be difficult if not impossible. 

Figures [12] and [13] show the spectral function (|85p for 
three different momenta: the M-point (face-center of the 
BZ) outside of the decay region, and the two points on 
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FIG. 12: (Color online) Magnon spectral function A(k, a;), 
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panel indicates position of the "bare" magnon peak. Solid 
arrows denote position of the renormalized quasiparticle peak, 
zero of Re[Gii(k,a;)~^], and the van Hove singularities. 



the TK line, k = (0.60857r, 0) and k = (0.27r,0), both 
inside of the decay region. The momentum for Fig. lisr a) 
is k*, which corresponds to the logarithmic peak in the 
on-shell Fk- The upper panels in Figs . [T^ and [T5I show the 
spectral function A{\i,uj), while the lower panels contain 
ReG"i(k,w), ImG'"i(k,w), and the two -magnon density 
of states (DoS): 

^k"''H=E'^('^-^1-^k-q). (87) 

q 

The energies of the non-interacting spin- waves ek are in- 
dicated by dashed arrows. The quasiparticle peaks in 
A(k, w) plots are marked by solid arrows. The lower 
panels shows that Re[Gii(k, w)~^] vanishes at the same 
UJ = Ek, as is expected for the pole behavior. For mo- 
menta inside the decay region, these peaks are also sig- 
nificantly broadened, see Figs. [THl^a), (b). While Ek differs 
very little from the one found in the numerical solution of 
the Dyson equation ([77)) . the damping in Figs.fT^a). (b) 
is somewhat smaller than the one found self-consistently 
from f77|) . As we discuss in Appendix [E] one can prop- 
erly define the quasiparticle residue even for the case of 
the non-zero damping. The issue of defining it at a sin- 
gularity point is also discussed there. 

Confirming our previous discussion, the spectral 
weight in the decay region. Fig. I13[ is non-zero above 
the two-magnon continuum oj-boundary. The lower pan- 
els in Figs. fTST a). (b) demonstrate that every singularity 
in the two-magnon DoS is reflected in both the real and 
the imaginary parts of the self-energy and, as a result, in 
the spectral function. For k = k* in Fig. fTST a) the bare 
magnon energy £k coincides with one such singularity, 
which corresponds to the decay into two magnons with 
energies e(k+G3)/2- This intersection with the singularity 
surface causes the anomaly in the on-shell spectrum, see 
Sec. IIVI For a non-singular momentum in Fig. [TWb'). 
which is not on the k* contour, similar singularity is 
above the magnon energy. 

Perhaps the most spectacular and also unexpected fea- 
ture of all the data in Figs.fTSTa). (b) are the sharp peaks 
at the bottom of the spectrum that are not associated 
with a concomitant peak in the two-magnon DoS. At 
the first glance, it may even be concluded that these 
peaks are the "true" , well-defined quasiparticle peaks 
with zero damping. A close inspection of Re[I](k, lu)] and 
Im[I](k, uj)] in the lower panels clearly connects the peaks 
in A(k, uj) to the log- and jump-like singularities in the 
one-loop self-energy. The origin of them is slightly more 
delicate than just a two-magnon DoS feature. The two- 
magnon DoS at the bottom of the continuum et in the 
decay region corresponds to the boundary to the emission 
of the £-|-Q magnon. Thus, with £f, = ek±Q, Aw = e — st, 
and q in the vicinity of ±Q, the threshold behavior of it 
is: 

dI"''"^{uj) (X J qdqSiAui-VQ\(i\)(xeiAuj)-{Auj) (88) 

Within the self-energy, the decay vertex exhibits anoma- 
lous behavior at small Sq — ±Q-fq for k away from k = 
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and from the decay boundary: f i(k, q) oc 1/ vRol- This 
yields a jump-like threshold behavior in the decay rate 



Im[E(k,cj)] cx J dqS{ALj-VQ\ci\) oc e(Aw), 



(89) 



and the concomitant log-singularity in Re[S(k, w)], as in 
Sec. |V1B. Thus, the weak threshold singularity of the 
type 0(Acl;) • Alu in the two-magnon DoS is enhanced 
by the singular decay vertex. This leads to a zero in 
Re[G'ii(k, o;)"^] and a "pseudo" -quasiparticle peak in 
A(k, w). We refer to these anomalies in Figs. [TWa). (b) 
as to the "edge" singularities. While in Fig. [TWa) this 
singularity is due to e_q magnon, in Fig. llSf b) there are 
two such singularities, one associated with £_q and the 
other e+Q magnon emission boundaries. 

As is discussed earlier, all such singularities should 
be regularized by the higher-order treatment. For the 
"edge" -singularities, the regularization requires summa- 
tion of the Pitaevskii's sequence, see Appendix [Cl since 
the external energy is purely real in A(\i,uj). 

To summarize, the conventional consideration of the 
single-particle spectral function in the triangular-lattice 
HAF is complicated by the logarithmic singularities asso- 
ciated with the van Hove singularities in the two-magnon 
continuum and with the singular behavior of the three- 
magnon coupling. Despite the contamination with these 
spurious features, a sensible analysis of the spectrum in 
terms of the broadened quasiparticle peak on the back- 
ground of the two-magnon continuum is still possible. 



VII. OTHER MODELS 

In this Section we demonstrate that anomalous fea- 
tures found in the spectrum of the triangular-lattice 
Heisenberg antiferromagnet are generic and appear in 
a wide variety of noncoUinear antiferromagnets. One 



straightforward generalization from the considered case 
is to the easy-plane XXZ model. Another modifica- 
tion consists of reducing lattice symmetry while keep- 
ing isotropic interaction between spins. This gives the 
so-called J-J' model on an orthorhombically distorted 
triangular lattice: the Heisenberg exchange along hori- 
zontal chains J is stronger than the interchain coupling 
J'. Another system discussed here is the kagome- lattice 
antiferromagnet with the XXZ anisotropy. In the follow- 
ing we outline basic kinematic conditions for two-magnon 
decays in these models and emphasize the singularities in 
the magnon spectrum. 



A. XXZ model 

The XXZ antiferromagnet on a triangular lattice is 
defined by the following Hamiltonian 



QX QX 



qV qV 



■aSfS' 



(90) 



For the easy plane system with the anisotropy parameter 
a = Jz/J < 1 spins form the same 120° structure as in 
the Heisenberg case. In the harmonic approximation the 
spin-wave energy is given by a simple modification of the 
Heisenberg formula (fTS]) : 



3J5V(l-7k)(l + 2a7k) 



(91) 



In addition, the three-boson interaction terms retain the 
same functional form as in Eqs. (fT9|l ~([22 |l . Therefore, 
for the 1/S consideration of the XXZ model the changes 
concern only Bogolyubov parameters and quartic terms. 
While a detailed consideration of this model is beyond the 
scope of the present work, we would like to highlight two 
ubiquitous features of its spin- wave spectrum determined 
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FIG. 14: (Color online) Magnon energy for the XY spin-1/2 triangular-lattice AF along representative directions. Dashed line 
is the linear spin- wave dispersion and solid line is including the first 1/S correction. 



by the three-boson interactions: strong renormahzations 
and decays. 

Let us first focus on the strongly anisotropic case and 
put a = 0. The ground-state energy for such an XY 
antiferromagnet is 



E^,,,/N = --JS' 



1 



0.064515 0.013326 



25* 



i2sy 



(92) 



in agreement with Ref . |33. The harmonic spin- wave spec- 
trum for this case is shown in Fig. [T3] by dashed lines. 
Because of the reduced spin-rotational symmetry there 
is only one acoustic branch near the F-point. The zero- 
point fluctuations reduce the sublattice magnetization to 
(S) = S — 0.051467 in the linear spin-wave approxima- 
tion. Even for the spin-1/2 case this amounts only to 
a 10% renormalization. Yet, the 1/5' corrections to the 
spectrum shown by solid lines in Fig. [TJ] are even larger 
than in the isotropic case! Magnon bandwidth narrowing 
in a = case is almost 50% of its bare value. One can 
also observe that the much-discussed "roton-like" min- 
imum at the M-point is much more pronounced here 
than in the spin- wave results for the Heisenberg limit. 
The origin of this minimum can be traced to the ID- 
like van Hove singularity in the two-magnon density of 
states at the Af-point, similarly to the isotropic case, 
Sec. lIVl B. In addition, one can notice that the top of the 
renormalized magnon band exhibits much more extended 
flat regions than in the Heisenberg limit as well as some 
weaker minima (between F and K) and other extrema. 
These features must affect the thermodynamics of this 
model substantially. Thus, interpolating between a ~ 
and a = 1 cases, one can conclude that the anharmonic 
three-boson terms lead to very strong spectrum renor- 
malization throughout the BZ for the XXZ model on 
the triangular lattice for all values of a. 

On the other hand, magnons in the a = XXZ 
model are stable at T = and decays are absent. To 
study the evolution of the two-magnon decays we begin 



with the nearly Heisenberg limit 1 
with q = ±Q are gapped with eq = 



a <C 1. Magnons 
This 



v3(l7«)/2 

has two immediate consequences in comparison with the 
a = 1 limit: (a) magnons at k — *■ cannot decay into 
two £q magnons up to k- values such that = 2eQ and 



a = 0.99 




a = 0.96 




FIG. 15: (Color online) Decay region and singularity lines for 
the XXZ triangular-lattice antiferromagnet with a = 0.99 
(top) and a — 0.96 (bottom). Definition of lines is the same 
as in Fig. [§] 
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FIG. 16: (Color online) Left panel: the 3D shape of the linear spin-wave frequency u)]^ = e-^/2JS for the J — J' antiferromagnet 
with J'/ J = 0.34. Central panel: the decay region and the singularity lines at J' / J = 0.34 (definitions are the same as in 
Fig. [9}. Right panel: same for J' / J = 0.7. 



(b) Q-magnons become stable themselves for the same 
reason. As a result, the star-shaped decay region of 
Fig. [9] develops a hole in the middle and has vertices 
shrunk and rounded, see Fig. [151 The evolution of the 
character of the decay boundary with a is non-trivial. 
Initially, the emission of a Q magnon remains an abso- 
lute minimum of the two-magnon continuum for most 
of the decay region except the tips of the star vertices. 
At ai « 0.993 the decay into a pair of non-equivalent 
magnons switches from being a line of saddle points into 
the absolute minima of the continuum and takes over 
the decay boundary, compare dotted and dashed lines in 
Figs. [H] and [TS] Without going into details, we simply 
note that in such a case the real part of the spectrum is 
singular, ReS(A;fc,cj) cx Inj/c — fcb|, see Eq. ([75)) . which is 
similar to the case of excitations in ^He near the thresh- 
old of decay into a pair of rotons.^^'**^ Similarly to the 
latter case, finding the spectrum near the boundary will 
require summation of the higher-order diagrams. Fig. [151 
shows the evolution of the decay region and singularity 
lines between a = 0.99 and a ~ 0.96. Further decrease 
of a completely eliminates the decay region at around 
a ~ 0.92. Therefore, magnon decays are present in an 
anisotropic triangular-lattice AFs, but only for not very 
strong anisotropics. 

B. J- J' model 

Another variety of generalizations of the nearest- 
neighbor Heisenberg model can be generated by spa- 
tial anisotropics, keeping spin-space intact. A particular 
model of this type corresponds to an orthorhombically 
distorted triangular lattice: 

X z.—z. 

n = J^S,-Sj+ J' Si ■ Sj , (93) 

(ij) {ij) 

in which interactions along the "ID chains" running par- 
allel to the X-axis is J, while zig-zag interaction between 



the chains is J'. This model has attracted a lot of atten- 
tion due to experimentally available systems, Cs2CuCl4^ 
and Cs2CuBr4''^ with J' /J is w 0.34 and « 0.7, respec- 
tively. For the former system, a comprehensive experi- 
mental neutron scattering analysis of the spin-excitation 
spectrum has been performed^ and an extensive theoret- 
ical analysis using both the spin-wave theory and the ID 
spinon-based approach have been carried out42ii^i^!^ 
We do not intend to repeat any of these calculations here, 
but would like to emphasize that a substantial broaden- 
ing of the spin-waves in a major part of the BZ must 
persist throughout the phase diagram of the J- J' model. 

The ground state of the classical J- J' model is an in- 
commensurate spin spiral. In the harmonic approxima- 
tion, the energy of the spin- waves in the J- J' model with 
J'/J < 2 isiii 

Ek = 6SJ^ (7k - 7q) ((7Q+k + 7Q-k) /2 - 7q), (94) 
where 

7k = - [cos kx + 2y cos — cos ^ j , (95) 

and the ordering vector Q = (Qx , 0) is given by 

Qx^Tr + 2 arcsin( J72 J) . (96) 

Since the 15*0(3) symmetry is preserved in the J- J' 
model, the Goldstone modes remain at k = and 
k = ±Q points with the incommensurate ordering wave- 
vector Q, Eq. ([96]) . Thus, because the velocities of these 
modes must generally stay different, magnon decays will 
always be allowed. In Fig. [T6| we show the 3D shape of 
the linear spin- wave energy for J'/J = 0.34. All nota- 
tions and symmetry points are the same as in Fig. [21 kx 
is in units of 1/a and ky is in units of 1/6, where a and h 
are the lattice constants in the chain and interchain direc- 
tions, respectively. One can observe a more sophisticated 
shape of the dispersion with three saddle points at dif- 
ferent energies (M, M', and between M' and K' points). 
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Comparing this Figure with the J = J' case of Fig. [2] 
one also finds that with the decrease of J' the spin-wave 
spectrum develops a low-energy branch in the y-direction 
(between the chains with strong J). That should make 
the rest of the spectrum prone to decays into it. 

To analyze the decays on a more qualitative level we 
show the decay regions and the singularity lines for two 
representative values of J'/ J in Fig. [111 First, the bound- 
ary of the decay region is determined by the emission of 
the gapless ±Q magnon for any value of J' / J as in the 
J' = J case. Second, the incommensurability of the or- 
dering wave-vector Q does not change the kinematics of 
the decays near k — > where the spin-waves can always 
decay into two modes near -|-Q and — Q. However, the 
incommensurability forbids the decays from near the Q- 
point into the vicinity of — Q point since the 2Q wave- 
vector is not equal to any reciprocal lattice vector any- 
more and quasi-momentum cannot be conserved in such 
a decay. This leads to shrinking of the vertices of the 
star-shaped decay region in the A:a;-direction, see Fig. [TOl 
However, the vertices expand considerably in the ky di- 
rection. The singularity lines due to the saddle points 
in the two-magnon continuum also expand and stretch 
in the ky direction as J' is lowered, in agreement with 
the expectations from the shape of the dispersion. Over- 
all, the decay region grows with the decrease of J' . At 
J' « 0.34J, relevant to Cs2CuCl4,^ the decay region cov- 
ers most of the BZ, see Fig. [TBI 

Within the spin-wave theory, a significant phase space 
volume for the decays found at J' / J = 0.34 can be some- 
what compensated by smaller noncollinearity of the spins 
which reduces the decay amplitudes. In real CS2CUCI4, 
additional Dzyaloshinskii-Moriya interactions make this 
angle large (close to 90°) between the spins in the nearest- 
neighbor ID chains,— but almost antiparallel within the 
individual chains with strong J . This means that, effec- 
tively, the decay vertices are proportional to the weaker 
coupling J'. However, the decays of the spin- waves still 
result in a substantial damping, see Refs. mill. Alto- 
gether, the J-J' model should exhibit magnon decays 
and singularities in their spin- wave spectrum throughout 
all the ranges of J' where the spin-wave theory is appli- 
cable. 



C. XX Z model on the kagome lattice 

The Heisenberg antiferromagnet on the kagome lattice 
is magnetically disordered at T = in both the classical 
{S = 00) and the 5* = 1/2 limits4^iiSi^ The degeneracy 
of the classical ground state reveals itself in the presence 
of the dispersionless zero-energy branch of magnons in 
the harmonic spectrum. The easy-plane anisotropy does 
not lift such a degeneracy, while opening up a constant 
gap for the zero-energy mode. A somewhat different 
behavior is realized in the recently discovered kagome- 
lattice compound potassium jarosite.^'^ In this antiferro- 
magnet, the Dzyaloshinskii-Moriya interactions lift the 



zero-energy mode to finite energies removing simultane- 
ously the classical degeneracy in favor of a so-called q = 
spin configuration. Still, the easy-plane XX Z model. 



■aSfS' 



(97) 



is more advantageous for a qualitative consideration of 
the magnon spectrum, because analytic expressions for 
the spin-wave energies can be easily derived. 

Since the unit cell of the kagome lattice consists of 
three atoms, there are three branches of magnetic ex- 
citations. For the 120° structure, the energies of these 
branches in the harmonic approximation are given by: 



where 



,(2,3) _ 



2^1 - a = luq , 



(98) 



(99) 



(I — a) ( I 

1 - a7k - 1 ± Vl + 87k 



with 



7k — cos k\ cos k2 cos fcs 



(100) 



and k\ — k^, ^2,3 = ±kx/2 + ky\/3/2, respectively. Thus, 
there is a gapped dispersionless mode, wq, gapped disper- 
sive one, uj'^K and gapless one, uj^\ see Fig. [T71 where 

w^*'* are shown for a = 0.95. In the real system, the 
lowest branch is weakly dispersive. ^'^ Using symmetry 
consideration, the BZ for the kagome lattice HAF can 
be reduced to the one smaller than the triangular-lattice 
BZ, see Ref. El In Fig. [13 TXY cuts are according 
to the notations of this work, with X = (0, 27r/3) and 
y = (7r/2,^/2V3). 

Since the three modes in (|99p are some linear combina- 
tions of local spin- flips, the three-boson interaction due 
to non-coUinearity of the spin structure will necessarily 
facilitate couplings of all three branches with each other. 
Such a coupling should not have any apparent smallness 
aside from the 1/VS factor, same as in Eq. (|20l) . There- 
fore, some qualitative conclusions about magnon decays 
in this system can be made without the detailed technical 
analysis. 

A distinct feature of the spin-waves in the kagome AFs 
is an (almost) dispersionless mode. For the decays this 
means that in the decays involving one of such magnons 
the momentum conservation can be ignored. This imme- 
diately implies that any spin-wave with the energy ex- 
ceeding ujo will have a finite life-time due to decays into 
one loq magnon and a Goldstone mode of the gapless 
branch. Thus, all the excitations above the horizontal 
dashed line in Fig. [T7] will be damped due to such decays. 

Another observation is that there must exist a particu- 

(2) 

larly strong singularity in both dispersive branches lo^^ 
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FIG. 17: (Color online) Linear spin-wave frequencies of the 
excitations in the XXZ model on the kagome lattice, Eq. (|99p 
for a = 0.95, along VXY cuts. Notations for the BZ are the 
same as in Ref. i49j. Lines and shaded area are the decay 
boundary (dotted), strong singularity line (solid), and the 
decay region, respectively, see text. 



and Ljj^ ' at the energy twice the energy of the dispersion- 
less mode. This is because the self-energy near 2wo will 
have a resonant-like form: Sk(w) cx: \y\^ j (uj — 2ujq). Even 
for the case when the lower branch has some residual dis- 

f2 3") 

persion, the damping of ' must be anomalously large 
near 2clIo in comparison with the rest of the spectrum. 
This should be valid even for large values of spin and 
in both 2D and 3D systems. Such a singularity is much 
stronger than the logarithmic singularities discussed in 
the rest of this paper and should be readily observed in 
experiment. 



D. Summary of Sec. IVlTl 

Summarizing the above examples, we conclude that 
strong renormalization, significant damping, and singu- 
larities must be common features of the excitation spec- 
tra of a wide variety of noncoUinear AFs. The enhance- 
ment of damping along certain contours due to singular- 
ities should be considered as fingerprints of magnon de- 
cays. This should help to distinguish the decay-induced 
spin- wave broadening from the other scenarios that yield 
broad spectra of the spin excitations. 



anomalous features of the spectrum that should be ob- 
servable in experiments are the substantial broadening 
of magnon peaks due to spontaneous decays in a large 
part of the Brillouin zone and strong deviations of the 
spectrum from the LSWT expectations for S — 1/2 sys- 
tems. The broadening should also be enhanced by the 
In S factor in the vicinity of certain contours in the mo- 
mentum space due to the van Hove singularities in the 
two-magnon continuum. Such an enhancement might be 
more visible in the systems with larger spins, such as 
the spin-5/2 triangular-lattice RbFe(Mo04)2,' although 
the overall damping will be smaller due to the 1/5' ef- 
fect. As the measurements of the lifetimes of spin excita- 
tions in the neutron-scattering experiments are expected 
to improve drastically in the future, this will allow for 
the detailed analysis of the damping. Therefore, the en- 
hancement of the damping along specific contours will be 
able to serve as a fingerprint of the spin- wave decays and 
could help to distinguish the spin- wave decay mechanism 
from the other scenarios such as fractionalization of spin 
excitations into spinonsi^i^ We have also demonstrated 
that the strong renormalization, decays, and singulari- 
ties in the spectrum are prominent in a wide variety of 
frustrated AFs of current interest. 

The strong cubic anharmonic coupling of the spin- 
waves induced by the non-collinearity is also the source 
of various technical challenges within the spin-wave ap- 
proach. Such challenges are largerly unfamiliar in the 
well-studied bipartite (coUinear) AFs where interactions 
are weaker and are of higher order. In this work, we 
have demonstrated that the spin- wave theory encounters 
some new problems in treating magnon interactions in 
the non-coUinear AFs. Generally speaking, the standard 
1/5 expansion for the spectrum becomes non-analytic, 
which is manifested by the singularities in the first-order 
1/5 self-energy. The origin of such singularities is related 
to the crossing of the single-particle branch with the sur- 
face of van Hove singularities in the two-magnon contin- 
uum. Regularization of such singularities also requires 
extra technical effort of going beyond the lowest-order 
1/5 approximation. 
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their anomalous spin- wave spectra. The highlights of the 



We are grateful to O. Starykh and A. Chubukov for 
illuminating discussions. We also thank H. Tsunetsugu 
for bringing Ref. ^ to our attention. Part of this work has 
been done at the Max Planck Institute for the Physics of 
Complex Systems, which we would like to thank for hos- 
pitality. This work was supported by DOE under grant 
DE-FG02-04ER46174 (A. L. C). The research at KITP 
was supported by the NSF under Grant No. PHY05- 
51164 (A. L. C.). 



25 



APPENDIX A: HARTREE-FOCK CORRECTIONS 

In the harmonic approximation the foUowing Hartree- 
Fock averages are nonzero for the triangular-lattice HAF 



k k 



7k 



7k 



hi 



A 



3^7^ 



{ataj) = 7kUkWk T XI 



E 
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UkWk 



7k 

4 V Wk 



(Al) 



These four constants can be expressed through combina- 
tions of three two-dimensional integrals 



ci 



Cl 



E 



(7k)' 

Wk 



Co = 1.5747334 



-0.1042539 , C2 = 0.3444458 



(A2) 



As is discussed in Sec. IIIII A.3. the quartic terms (p4)) 
yield a correction to the ground-state energy that is given 
by the 4-boson averages. In the leading order, they can 
be decoupled into the bilinear combinations (|Aip using 
the Wick's theorem. The corresponding terms in Eq. 
are given by: 



i^ + m^ + A"^ , (A3) 
{a\aiaiaj) — 2nA -|- m5, {a^jOjajai) = 2nm + AS. 



This yields the ground-state energy correction from the 
quartic terms: 



A<5 



(A4) 



After some further algebra the above expression is con- 
verted into a combination of c; constants given by 
Eq. ([26]). 

A similar mean-field decoupling procedure yields the 
following correction to the harmonic part of the Hamil- 
tonian in terms of the original bosons: 



Mk = --J[2(co+ci-2c2-l)-h7k(co + C2-l + ici)], 



5Sk = ^ J [7k (1 + C2 



Co 



(A5) 



After that one uses the Bogolyubov transformation ^ 
and obtains the following quadratic form of magnon op- 
erators: 



(A6) 



The coefficients in this expression are related to 6Ay^ 
and ^i?k by 

£k ' = ("k + vl) ^^k - 2wkWk^Bk , 



B 



k 

(4) 



(ul + vl) (5Bk - 2ui,v^SAk , (A7) 



which leads finally to Eqs. 



and 



APPENDIX B: SUBLATTICE MAGNETIZATION 

Calculation of the second-order correction to the sub- 
lattice magnetization (|46p requires evaluation of the 
lowest-order contributions to {bl^bk) and (6k6-k), which 
are expressed in terms of the normal and the anomalous 
Green's functions, respectively: 

{bibi.)=i I ^Gn(k,c.)e-^ , 

^Gi2(k,c^). (Bl) 

One needs to keep only the first-order terms in the pcr- 
turbative expansion of the Green's functions. The di- 
agonal average (^k^k) has a single nonzero contribution 
determined by ^Gii(k, cj) = G§(k, a;)S^'^^ (k, w): 

which yields 

+ H W)' ,33) 



where we have transformed to the dimcnsionless vertex 
and frequencies in the last expression. 

The off-diagonal average (6k&-k) is determined by the 
lowest-order anomalous self-energies: 

SGuiKu;) = Go(k,c^)Go(-k,-o.)I]i2(k,o.) (B4) 

with Si2(k,cj) = I]HF(k) + si2Hk,tj) + Y.[f(k,uj). 
The first frequency-independent contribution from the 
Hartree-Fock self-energy ([57)1 yields 



7k(l - 7k) /I , , 



9 9 , 

= -325'^'^+32^('^"'^)E 



7k(l - 7k) 



Two other terms give identical contributions to (6k^-k) 
with the net result 



ss" = A V ^ V ri(q;k)r2(q;k) 



(B6) 
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Combining all of the above terms together, we obtain: 



SS2 = -^C1C2 + ^(C2-C1)^ 

k 



7k(l - 7k) 
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9 7k 

k k q 



ri(k,q)r2(-k,q) 



- t^k-q + t^k 



(B7) 



1 



27k 



r2(k,q)^ 



as a double integral over k and q using Eq. (|B8[) . An ad- 
ditional technical detail is that one should use projector- 
type multipliers Pi — |(1 — 7k) and P2 = \{1 + 27k) 
to avoid introducing extra singularities by the above 
conversion in the second term. These projectors obey 
Pi + P2 = 1 and guarantee convergence of the integrals 
close to the k = and k = Q points, respectively. Alto- 
gether, this has led us to the analytically identical, but 
numerically regular form of the l/S*^ correction to the 
on-site magnetization: 



As discussed in Sec. IIIII B.2. one cannot use Eq. (jB7[) 
for numerical evaluation of 6S2 directly as it leads to 
ambiguous results or simply does not converge. The way 
to regularize the above expression is to use an analytical 
insight-^"'^ The integrand in the third term in (|B7p can 
be reduced precisely to the divergent second term for 
k^ Q: 



9 7k ^ri(k,q)r2(-k,q) 



4 uol ^ 

K q 



- ^k-q + ^k 



= -Y^(c2-ci) 



k^Q 

7k(l - 7k 



+ OK^). (B8) 



Then, the proper subtraction of the leading singularities 
is ensured by re-expressing the second term in Eq. (IB7|) 



5^2 



16 



C1C2 -I- 5S 



3,1 



3,2 



8S- 



3,3 



(B9) 



where the second and third terms are regular from the 
start due to projectors and the last term is the regularized 
combination of all singular terms. Specifically: 



(1 + hk) 

Wk 



(| + 7k)E 



r2(k,q)^ 



(Wq -I- Wk+q + Wk)^ 



^ ^ CJq + t^k-q + t^k 

The combination of the divergent terms is given by: 



^^3,2 ^ o 7k 



J 



8S- 



3,3 



7kj 



E 

q 



7k 



ri(k,k-q)r2(-k,q) ri(Q,Q-q)r2(-Q,q) 



t^q + Wk-q + Wk 



2c^k(l + hk)r2(k,q)'' 

3 (tJq + Wk+q + t^k)^ 



The results of the numerical integration of the individual 
terms are: 

-— C1C2 = 0.20199272 , 5S^'^ = 0.017918(1) , 
16 

55^3^2 ^ 0.025496(2) , 5S^'^ = -0.074660(5) . 

Altogether, they lead to the following value of the second- 
order correction 

5S2 = -0.011045(5) . (BIO) 

APPENDIX C: HIGHER-ORDER 
SINGULARITIES 

Singularities in the perturbative calculations of the 
bosonic excitation spectra are known since the early 
works by Pitaevski i^^i"^^ on the termination point in the 
phonon branch of ^He. There is a renewed interest in the 
similar problems in the context of various spin systemSf22. 
On the other hand, the logarithmic singularities in the 
fermionic spectra are also known to occur since the early 



works on the edge-singularities in metals;^ in the context 
of numerous aspects of ID Luttinger liquids)^ as well as 
in the newer physical systems such as graphenci^ 

In all these problems, some physical process, often of a 
threshold nature, leads to a non-analytic behavior of var- 
ious quantities. Such a non-analyticity manifests itself in 
a breakdown of the perturbative expansion, i. e. as a sin- 
gularity. The theoretical challenge is to reconstruct the 
original non-analytic behavior from the singular terms 
in the perturbative expansion. In the Pitaevskii's case, 
such a reconstruction is rather straightforward and con- 
sists of resummation of the leading divergent terms of the 
"ladder" (RPA) type. In such a case, the original non- 
analyticity is straightforwardly related to the singularity, 
e.g. l/ln|(jj| to ln|ti;|, respectively. In the fermionic sys- 
tems, such a resummation is more complicated and in- 
volves the infamous "parquet" diagrams. However, often 
enough, the reconstruction of the original non-analytic 
behavior is still possible from the analysis of a few most 
divergent terms of the expansion. In some well-known 
cases, the log-singularity resumes in a non-analyticity of 
a non-trivial fractional power-law type, |u;|"i^ 
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(a) 

(b) 

(c) 
(d) 





FIG. 18: (Color online) Schematic picture of all topologically 
inequivalent diagrams of order 1/S, (a), and , (b)-(d), re- 
spectively. The leading divergence in l/S diagrams is In | Aa;| 
(first diagram in (a)), in order it is (ln|A[ij|)'^ (all dia- 

grams in (d)). First two diagrams in (d) are identical. 



Since in our case we deal with bosonic excitations, one 
can expect that the Pitaevskii's consideration is the most 
relevant one. While this expectation, with some minor 
corrections, turns out to be true, we would like to outline 
some differences and similarities of the singularities in 
the spectra of noncoUinear AFs in a somewhat broader 
context. 

One obvious difference of our problem from the 
fermionic case is the presence of two coupling constants, 
three-particle, V3, and four-particle, V4, depicted in 
Fig. [3] Because of that, the diagrammatic expansion 
is more complicated in our case. Since the singularity 
occurs for any value of the spin, it is natural to group 
the diagrams by their order in 1/5. By construction, 
the three-particle vertex is of order 1 / ^/S relative to the 
magnon energy £k and it must occur in pairs in any self- 
energy diagram, while the four-particle vertex is of order 
1/S. We depict schematically all topologically different 
diagrams of order 1/5 and 1/5^ that occur in such a 
theory in Fig. [TSTa) and Fig. [T^b)-(d1. respectively. 

The singularity in our case is due to the "bubble" di- 
agram [#1 in Fig. [T5f a)]. which leads to the logarithmic 
term for certain k- values of external lines of the diagram. 
A simple analysis shows that out of the 1/5^-terms only 
diagrams in Fig. [TST d) yield the higher power of the log, 
(In |Aa;|)^. In fact, the first two diagrams in the (d) group 
are identical and can be seen as vertex corrections to the 
bubble in Fig. [TSTa). The last diagram in the group (d) 
is the standard "two-bubble" member of the ladder se- 
quence considered by Pitaevskii. A closer inspection of 



the second diagram in the (d) group in Fig. [18] makes it 
clear that it is also the "two-bubble" diagram with the 
virtual magnon line and two 3-magnon vertices joining 
the "bubbles" . This modifies consideration by Pitaevskii, 
albeit only quantitatively, as it introduces extra interac- 
tion together with some retardation into the "ladder" 
sequence. To rule out the presence of any "non- ladder" - 
type diagrams with the higher powers of the log and, 
hence, to eliminate the possibility of more complicated 
form of the non-analytic behavior, the following observa- 
tion is useful. In the parquet consideration^^ it is both 
particle-particle and particle-hole diagrams that are sin- 
gular. In our case of bosons with non-conserved number 
of particles at T = 0, the particle-hole diagrams are all 
identically zero. Thus, the leading power of the log of any 
given diagram is simply the number of non-equivalent re- 
gions of the diagram in which a vertical cut intersects 
only two magnon lines. Therefore, the most-divergent 
diagram sequence is comprised of the "ladder" -like dia- 
grams only, with both V4 and (Va)^ between the "bub- 
bles". This can be confirmed by an explicit analysis of 
the "dangerous" diagrams in the next, 1/5^ order (not 
shown in Fig. [T5)l . 

Another qualitative difference of our problem from the 
fermionic case in which the "parquet" summation is suc- 
cessful is that the latter consideration often exploits the 
long- wavelength character of the problem to obtain a uni- 
versal answer. In our case, the singularity occurs at short 
wavelength. Therefore, even in the hypothetical gener- 
alization of our problem to the bosons with conserved 
number of particles and singular particle-hole "bubbles" , 
the "parquet" -like analysis is unlikely to be successful as 
the coefficients of the leading log-powers will be indepen- 
dent from each other. 

Altogether, presented analysis of higher-order diver- 
gences in the diagrammatic sequence for the magnon 
propagator in the noncoUinear AFs shows that the sin- 
gularities are regularized by summation of the "ladder" 
diagrams of the Pitaevskii type with the minor modifica- 
tion of the interaction vertex. The regularized result can 
be written schematically as: 



E cx 



V^U , B 



1 - vn 



(Cl) 



where V = V4 + V3GV3, and 11 is the singular "bubble" 
contribution (^ In |Aa;|). 



APPENDIX D: OFF-SHELL DYSON EQUATION 
FOR COMPLEX ENERGIES 



Suppose that the Dyson equation 

w-Wk-X;(k,w) = 0, (Dl) 
has a solution at a complex energy 

a; = (5u - iFk, Fk > , (D2) 
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where the sign of the imaginary part is dictated by causal- 
ity. The Dyson equation can be written for the real and 
imaginary parts as 

iik = ^k + Rc[I](k,(ik-irk)], (D3) 
Tk = -Im[S(k, - tTk)] > . (D4) 

In the standard textbook approach^ the damping is as- 
sumed to be small, Fk ^ Wk, and is neglected from the 
self-energy in the right-hand side of Eqs. (|D3p and (|D4p . 
In that case, only the real part of the Dyson equation 
needs to be solved self-consistently while the damping 
is simply evaluated from the imaginary part of the self- 
energy at the real w = Wk, which is the solution of (|D3|1 . 

Since we have encountered a problem in which the 
imaginary part of the self-energy diverges at some k- 
values for lo along the real axis, we would like to de- 
viate from the standard approach and keep finite Fk 
in both parts of the Dyson equation, solving them self- 
consistently for both Wk and Fk- However, on this path 
one discovers a difficulty that is unrelated to the singular- 
ity and is much more generic. We believe that it deserves 
a separate discussion. 

Let us consider the self-energy of very general form: 



S(k,c.)=^. 



Ek-c 



iO 



(D5) 



The self-energy of such type appears in the problem 
of phonon interaction with electron-hole continuum, 
Cherenkov radiation, pair production in QED, to name 
a few. All of these problems involve a (virtual) decay 
of a particle into two other particles. For simplicity, the 
products of the decay are assumed to be free particles (no 
damping for eq's). The difficulty arises when substituting 
w = Wk — iFk in the denominator of (|D5[) . This shifts the 
pole in (jPip into the wrong (upper) w-half-plane making 
Fk < 0, which violates causality and renders the Dyson 
equation unsolvable. 

The resolution of this problem within the formalism 
of causal diagram technique requires one to define prop- 
erly the energy conservation for particles with finite life- 
time. The self-energy (jDSp is the result of integrating 
the "bubble" -like diagram over the internal frequency. 
Specifically, assuming for simplicity that the decay is into 
particles of the same species, one can write: 



l^k,qP 



(w-e-Ek-q+iO) {e-Sq+iO) 



(D6) 



where = J de/2TT. This step seems to make the prob- 
lem even worse: substituting cj = Wk — iFk in (jD6p shifts 
the pole of the first Green's function into the wrong half- 
plane, which renders the integral over internal frequency 
zero. 

Making one more step back we recall that Eq. (|D6|) 
comes from the standard diagrammatic expression that 



explicitly reflects the conservation of energy within the 
decay process: 



1 (ei-£k-q+«0) {e-Eq+iO) 



(D7) 



Integrating (jD7P over ei gives (|D6p . Now, the only place 
where "external" energy uj enters the self-energy is the S- 
function. Thus, one needs to generalize the 5-function in 
Eq. (jD7P for the case when u; turns complex. An obvious 
trick is to represent the ^-function by the Lorentzian of 
width Fk. Then one can rewrite S(k, lo) for uj = Wk — *Fk 
as: 



S(k,L^) 



|Vk.qP2Fk 

Jiik-ei-£)2+F2 
1 



r- (E)8) 

(ei-£k-<i-l-iO)(£-eq+iO) 

Note, that the above expression does respect the causal- 
ity. Integrating it by standard means we obtain the de- 
sired well-behaved result: 



I](k,c.)=^^ 



^k - £q - Ek-q + «Fk 



(D9) 



One can see that, starting from the basic diagrammatic 
rules we obtain the result for the self-energy with the 
"wrong" sign of the imaginary part of the energy of the 
decaying particle. This resolves the problem. In the form 
(ID9p the self-energy has the quasiparticle pole in the cor- 
rect (lower) half-plane and yields the correct sign of Fk in 
the Dyson equation. Thus, the prescription is: the energy 
of the decaying particle in the self-energy (|D5p should be 
taken from the advanced Green's function. The "proper" 
Dyson equation then reads as 



LJ — LJk — S(k, UJ*) — 0, 



(DIO) 



for complex lu = uj^ — i^k- 

Using the Matsubara technique, the problem is re- 
solved without any tricks with the (5-function, but sim- 
ply by forcing all the poles to be in the correct half-plane. 
This does not give any additional insight, but means that 
the retarded self-energy simply corresponds to the com- 
plex conjugation of w: Sret(k, Wk — iFk) = S(k, cjk+iFk). 
Thus, one can write the Dyson equation in an explicitly 
self-consistent form as: 

t^k — «Fk — t^k — S(k, cjk + «Fk) = (F>11) 

with I](k, cij) from (|D5p . Changing Fk to an infinitesi- 
mal S, this expression is nothing but the standard Dyson 
equation obtained from the Matsubara approach, see 
Ref.lH, Eqs. (3.140)-(3.142). 



APPENDIX E: QUASIPARTICLE RESIDUE 

Here we discuss an extension of the definition of the 
quasiparticle residue to the case of a particle with a fi- 
nite lifetime. The concept of the quasiparticle residue 
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is introduced in the context of the problem of interact- 
ing particles where it is assumed that after "dressing" the 
Green's function retains a well-defined pole at low energy. 
Close to that pole one can write the Green's function in 
the main-pole approximation: 



Gk(^) 



G 



'incoh 



(El) 



where the incoherent part of the Green's function 
Q^ncoh^^^ is regular at the "new" quasiparticle energy 
w = cDk and Zk is the quasiparticle residue, Zk < 1. 
is also the weight associated with the (5-functional peak of 
the Green's function at lu = uj^. We wish to extend this 
to the case when the quasiparticle has a finite lifetime. 
In other words, we hope to be able to write: 



Gk(t^) 



1 



- Wk - Sk(^^) 



(E2) 



Wk -I- iTk 



incoh 



To get from the left-hand side to the right-hand side we 
need to assume that there exist a solution of the Dyson 
equation 



w — o^k — Sk(ti') = 

at some complex a; = cjk ^ *rk, where 

^^k = ^k + Re[Ek(wk - «rk) 
Tk = -Im[i;k(i)k - irk)] • 



(E3) 



(E4) 
(E5) 



Then, as in the standard approach, we proceed by adding 
and subtracting Sk(t3k — irk) = Wk — i^k — i^\L to the 
denominator of Gk('^): 



Gk(c^) 



1 



- t^k - Sk(ti') 



(E6) 



1 



w - cZik + jFk - (Sk(w) - Sk(a;k ~ iFk)) 
This can be rewritten without approximation as: 
1 



Gk(u;) 



Wk + «rk 

Sk(t^) - 5]k(wk - irk) 



(E7) 



Near the pole, 
yields: 



CJk 



Gk(a;) 



cj — cjk + iFk 
iFk, this expression finally 

(E8) 



where 



Zv = 



1 - 



— tJk + i^v 



ask(co) 



(E9) 



which coincides with the standard definition of Zk up 
to the change of Wk to Wk — irk.^** Therefore, generally 



speaking, the quasiparticle residue obtained this way is 
complex. However, one would prefer to have Zk real 
in accord with the expectation that the area under the 
Lorentzian (broadened (5-peak) should give the quasipar- 
ticle weight. For that, let us consider the spectral func- 
tion: 



^k(w) 



1 / rkRe[Zk] 
{ll> - Wk)Im[Zk' 



■(ik)2 



(ElO) 



where the first term has the Lorentzian form while the 
second one vanishes at Wk = Wk- Clearly, the second part 
of the expression is odd in w and the spectral weight of 
the quasiparticle peak, which is given by the integral of 
A]s,{uj) along the real axis: 



Ak(w)dw = Re[Zk] 



(Ell) 



is simply the real part of Zk in (jE9p . This gives a proper 
definition of the "generalized" quasiparticle residue: 



Zk = Rc 



1 - 



aSk(c^) 



duj 



(E12) 



Interestingly, this definition differs from the one in 
Ref. M- 

Our interest in this problem is twofold. First, we deal 
with quasiparticles that have finite damping due to de- 
cays. Second, the decay part of the self-energy, Eq. 
near the log-singularity has derivatives that are even 
more singular. In particular, an attempt to calculate the 
quasiparticle residue using the on-shell approach near the 
saddle-point singularity k ^ k* leads to: 



duj 



1 



ri(q;k)2 



2 ^ (cJk - 



■ Wk- 



whose imaginary part is divergent as oc (fc — fc*)~^ and 
the real part has a (5-function singularity at fc = fc*. If, 
on the other hand, one uses the solution of the Dyson 
equation together with our more general definition of Z^ 
(jE12|) . the result becomes regular. This is yet another 
way of saying that the l/S'-expansion in the noncoUinear 
AFs is singular and the usual on-shell approach cannot 
be used. 

Interestingly, the off-shell consideration gives, after 
some algebra, that in the large-S* limit, at the singu- 
lar k-points {dT,/duj) cx l/ln5. This means that the 
quasiparticle residue reaches the classical limit (=1) very 
slowly as 



Zw 



A 



at these points. 
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